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Abstract: 

We calculate the complete 0(a) electroweak radiative corrections to e + e~ — > WW — > 
4/ in the electroweak Standard Model in the double-pole approximation. We give analyt- 
ical results for the non-factorizable virtual corrections and express the factorizable virtual 
corrections in terms of the known corrections to on-shell W-pair production and W de- 
cay. The calculation of the bremsstrahlung corrections, i.e. the processes e + e~ — > 4/7 in 
lowest order, is based on the full matrix elements. The matching of soft and collinear sin- 
gularities between virtual and real corrections is done alternatively in two different ways, 
namely by using a subtraction method and by applying phase-space slicing. The O(o) 
corrections as well as higher-order initial-state photon radiation are implemented in the 
Monte Carlo generator RacoonWW. Numerical results of this program are presented 
for the W-pair-production cross section, angular and W-invariant-mass distributions at 
LEP2. We also discuss the intrinsic theoretical uncertainty of our approach. 



June 2000 



1 Introduction 



At present, the focus of electroweak Standard Model tests lies on W-boson-pair pro- 
duction at LEP2, e + e~ — > WW — > 4/ [ fj, @, 0]- While the invariant-mass distributions 
of the final-state fermion pairs allow for precise measurements of the W-boson mass, the 
total cross section and the angular distributions can be used to obtain information on 
the non-abelian couplings of the W bosons. LEP2 provides us with (9(10 4 ) W-boson 
pairs, thus leading to a typical experimental accuracy of 1% [ BJ. At a future e + e~ linear 
collider, the W-boson yield could be increased by two orders of magnitude [ |J], and the 
experimental accuracy will reach the level of some 0.1%. 

The Monte Carlo generators used in the past [ |5|] for W-pair production are based 
on the lowest-order matrix elements and typically include only the universal electroweak 
corrections. These comprise the running of the electromagnetic coupling, corrections as- 
sociated with the p parameter, the Coulomb singularity close to threshold, and photonic 
corrections in leading-logarithmic approximation. For the total cross section these gener- 
ators have a precision at the level of 1-2% at LEP2 energies [ || ^ §]; for distributions 
the accuracy in general is worse. While this was sufficient for the analysis of LEP2 data 
in the past, it becomes already insufficient for the present analysis and will be definitely 
inadequate for the final LEP2 analysis. Since the non-universal corrections increase with 
energy, the accuracy of generators including only universal corrections is at the level of 
5-20% [ H in the total cross section and thus not adequate at linear collider energies. 

In this paper we aim at a theoretical accuracy of < 0.5% at not too large energies^, 
which is sufficient for LEP2 and reasonable for a 500 GeV linear collider. This requires, 
on the one hand, to include the complete set of lowest-order diagrams for e + e~ — > 4f, 
which is already implemented in many existing generators. Moreover, the complete 0(a) 
corrections for e + e~ — > WW — > 4/ have to be taken into account, and care must be taken 
that also the leading higher-order corrections are included. 

The full treatment of the processes e + e~ — > 4/ at the one-loop level is of enormous 
complexity. Nevertheless, there is ongoing work in this direction [ [|. Moreover, the 
requirement to include the finite width of the W bosons poses severe theoretical problems 
with gauge invar iance. 

An economic approach for the calculation of the 0(a) corrections to e + e~ — > WW — > 
4/ consists in using the double-pole approximation (DPA). In DPA, only those terms are 
kept that are enhanced by two resonant W propagators and thus by a factor Mw/Pw 
with respect to all other contributions. With respect to the leading contributions, the 
corrections to the non-doubly-resonant contributions are typically of the order a/ir x 
Tw/^fw- Even when taking into account a conservative safety factor, this approximation 
is sufficient for our precision tag of 0.5% as long as the doubly-resonant contributions 
dominate the cross section. This is the case at LEP2 energies sufficiently above threshold 
if forward-scattered electrons or positrons in the final state (single-W production^) and, 
for processes involving the neutral current, fermion-antifermion pairs with low invariant 
masses are excluded. At higher energies, where diagrams without two resonant W bosons 

1 Above 0.5-1 TeV at least the leading electroweak logarithms at the two-loop level should be taken 
into account. 

2 A discussion of state-of-the-art calculations for single-W production can be found in Ref. [§. 
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become more sizeable, it may be necessary to suppress those by appropriate cuts, e.g. 
invariant-mass cuts, on the final-state fermions. These types of cuts must be applied 
anyhow in order to extract the physically interesting WW signal. 

The DPA has several advantages as far as the virtual corrections are concerned. Be- 
sides reducing the number of diagrams drastically, it naturally provides a gauge-invariant 
answer. In contrast to the corrections to the full process e + e~ — > 4f, the virtual cor- 
rections in DPA depend on the fermion flavours of the final state in a universal way, 
i.e. there is only one generic set of diagrams for all final states. Moreover, the invariant 
functions that build up the corrections do not depend on the decay angles but only on 
the production angle of the W bosons, thus allowing to speed up their calculation in the 
Monte Carlo generators considerably. Finally, the virtual factorizable corrections in DPA 
can be directly obtained from the existing results for on-shell W-pair production [ [10, |11 



and W-boson decay [ 0, |13 . The virtual non-factorizable corrections, i.e. those virtual 



corrections that cannot be associated with either W-pair production or W-boson decay, 



have a simple structure [ [14], [15|, [16], [17 



On the other hand, the definition of the DPA is problematic for the bremsstrahlung 
corrections. Photon radiation from the initial and from the final state leads to W propa- 
gators that become resonant at different locations in phase space. The situation depends 
on the regime of the photon energy. If the radiated photons are hard (E 1 ^> Tw), these 
locations are well separated, and the photon radiation can be unambiguously assigned 
to the W-pair production subprocess (at least theoretically) or to one of the W-decay 
subprocesses. For soft photons (E 7 <C Tw), the resonances coincide, and the DPA is 
identical to the one without photon. However, for semi-soft photons (E 1 ~ Tw), the 
resonant propagators overlap, and it is not obvious how the DPA can be applied without 
omitting or double-counting doubly-resonant contributions. For these reasons, we have 
preferred to employ the exact e + e~ — > 4/7 matrix element for the real bremsstrahlung 
corrections and the exact e + e~ — > 4/ matrix element for the soft and collinear singular 
virtual corrections, and to apply the DPA only to the non-leading part of the virtual 
corrections. Furthermore, all observables for real photons are based on a full lowest-order 
calculation for e + e~ — > 4/7 []. 

A possible strategy for a DPA for the 0(a) corrections to pair production of unstable 
particles was already proposed in Ref. [ |IBf. Recently different versions of DPAs have 
been used in the literature. 

A first complete calculation of the O(a) corrections for off-shell W-pair production, 
including a numerical study of leptonic final states, was presented in Ref. [ [19| . In that 
work the DPA is applied both to the virtual and real corrections using a semi-analytical 
approach. Moreover, also the phase space is treated in DPA. The proposed DPA for 
the bremsstrahlung process is based on the distinction between hard, semi-soft, and soft 
photons, without a common treatment of photons of these domains for all observables. A 
realistic photon recombination for collinear photon emission is not considered there. For 
instance, the invariant masses of the W bosons are reconstructed from the momenta of the 
corresponding decay fermions alone, without recombining soft and collinear photons. The 
resulting large shifts in the peak position of the W line shape, namely —20 MeV, —39 MeV, 
and —77 MeV for r + u T , and e + v e final states at a centre-of-mass (CM) energy of 



3 A survey of state-of-the-art calculations for e + e — > 4/7 is also contained in Ref. 
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184 GeV, respectively, result from mass-singular logarithms of the form a ln(m;/Mw) and 
are due to the absence of collinear photons in the definition of the invariant masses. 

0(a) corrections for W-pair production in DPA have also been implemented in the 
Monte Carlo generator YFSWW3 [ [SJ |]|. In Ref. [ |0j the 0(a) corrections to on-shell 



W-pair production have been combined with the exponentiation of photonic corrections 



from the initial and intermediate (WW) states. In Ref. [21], the final-state corrections 
have been added in the leading-logarithmic approximation. The normalization of the 
W decays is fixed using the corrected branching ratios. The non-factorizable corrections 
are not fully included but only as an approximation in terms of the screened Coulomb 
ansatz [ |22| . The spin correlations between the two W decays have been neglected in the 



contributions of the non-leading corrections. In Ref. [ the results of Ref. [ have 
been confirmed qualitatively. Moreover, it has been found that the shifts of the peaks in 
the W-invariant-mass distributions are reduced by taking into account collinear photons 
in the definition of the invariant masses, which effectively replaces the mass-singular 
logarithms by logarithms of a minimum opening angle for collinear photon emission. 

Another approximate version of the DPA was presented in Ref. [ ^3] and discussed 
for a purely hadronic final state. Both virtual and real corrections were treated within 
DPA, but non-factorizable corrections were neglected. The main emphasis was put on 
linear-collider energies, in particular on the comparison to a high-energy approximation [ 
24 1 for the one-loop correction. 

In this paper we present the first complete calculation of the 0(a) corrections for 
off-shell W-pair production in the DPA that has been implemented in a Monte Carlo 
generator. Numerical results of this generator, which is called RacoonWW, have already 



been presented in Refs. [ 25, pfl and [ 26] for LEP2 and linear-collider energies, respectively. 



The lowest-order cross section is calculated including the complete diagrams for any four- 
fermion final state or in the CC03 approximation, i.e. including only the doubly-resonant 
W-pair production diagrams. The complete virtual corrections to W-pair production and 
W decay and the virtual non-factorizable corrections are included in the DPA. As far as 
the real corrections are concerned, the matrix elements for the minimal gauge-invariant 
subset of e + e~ — > 4/7 including the doubly-resonant contributions corresponding to W- 
pair production (the CC11 subset) are included. The real corrections and the virtual 
corrections are matched in such a way that soft and collinear singularities cancel as far as 
they should. The treatment of these singularities has been implemented in two different 
ways, one of which uses the subtraction method described in Refs. [ ^7], |T7[] , the other one 
uses phase-space slicing. We use the exact four-fermion phase space throughout. Finally, 
we note that all parts of the calculations have been performed in two independent ways. 

The paper is outlined as follows: in Section Q we describe the general strategy of 
our calculation. In Section ^| details on the virtual factorizable and non-factorizable cor- 
rections in DPA are given. In Section |] we explain the matching of soft and collinear 
singularities between virtual and real corrections, both for the subtraction method and 
for phase-space slicing. The inclusion of higher-order initial-state radiation is described 
in Section [5|, and Section ^ contains a discussion of QCD corrections. Numerical results 
are presented in Section |7|, including predictions for the W-pair production cross section, 
angular and W invariant-mass distributions at LEP2, a detailed discussion of the intrin- 
sic ambiguities of the proposed DPA, and comparisons to results of other authors. Our 
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conclusions are presented in Section Q and the appendices provide some useful explicit 
results. 



2 Strategy of the calculation 

We consider the process 

e + (p + , ( x + )+e-(p_,a_) - W+(A; + , A+) + W"(A;_, A_) 

-> fx(h^x) + Mh,a 2 ) + f 3 (k 3 ,a 3 ) + h(h,a A ). (2.1) 

The arguments label the momenta p±, hi and helicities cr,j = ±1/2, Aj = 0, ±1 of the 
corresponding particles. We often use only the signs to denote the helicities. The fermion 
masses are neglected whenever possible. As a consequence, the helicities of the incoming 
electrons and positrons are related as a = o~_ = — a + (in the absence of collinear photon 
emission). For the W-pair-mediated diagrams, and more generally for the graphs of the 
CC11 class, the helicities of the outgoing fermions are fixed, a^ 3 = — cr 2 ,4 = —1/2, owing 
to the left-handed coupling of the W bosons. In general this does not hold for other 
background diagrams. 

The lowest-order cross section is calculated using the complete lowest-order matrix 

elements A^Bom^ 4 ^ °f R- e ^ \ 1^8 • 



The radiative corrections to (2.1) consist of virtual corrections, resulting from loop 
diagrams, as well as of real corrections, originating from the process 

e + (p + ,<7 + ) + e-(p_,<7_) - W + (A; + ,A + )+W-(^,A_) ( + 7 ) 

-> fi(h, vi) + /2(Afc, o"2) + f 3 (h, 0-3) + U(ki, 04) + l{k, A 7 ). 

(2.2) 

Both have to be combined properly in order to ensure necessary cancellations of soft and 
collinear singularities. 



2.1 Virtual corrections 

We treat the non- leading virtual corrections in double-pole approximation (DPA), i.e. 
we take only those terms into account that are enhanced by two resonant W propagators 
(doubly- resonant corrections). 



2.1.1 Virtual factorizable corrections 

In the DPA, there are two types of contributions, the factorizable and the non- 
factorizable ones. The former are the ones that can be associated to one of the pro- 
duction or decay subprocesses, the latter are the ones that connect these subprocesses. 
This splitting cannot be completely done on the basis of Feynman diagrams; otherwise 
the splitting would not be gauge-invariant. To define this splitting, we start by inspect- 
ing only those graphs that retain two resonant W propagators outside the loops, i.e. all 
diagrams of the generic structure shown in Figure [p. Since the loop corrections in such 
graphs can be associated either to the production of the W-boson pair or to the decay of 
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On-shell production On-shell decays 

Figure 1: Diagrammatic structure of virtual factorizable corrections to e + e" 



WW-+4/ 



one of the W bosons, this subset of graphs can be associated with the virtual factorizable 
corrections. The corresponding amplitude is of the form 



M 



R{k 



2 p 
_|_, 



(4-Af& + iM w r w )(A2 



M%j + i.lAvh 



(2.3) 



w 1 wj 



where k\ are the invariant masses of the virtual W bosons, 9 represents all other kine- 
matical variables, M\y is the renormalized W-boson mass, for which we take the on-shell 
mass, and Tw is the width of the W boson. 

The diagrams of Figure [l] do not form a gauge- invariant subset. However, in DPA 
a gauge-invariant contribution can be extracted by replacing the numerator with the 



gauge-invariant residue [ |29 



M 



DPA 



+WW->4/ 



U'i + iM 



w r w )(fc_ 



M$f + iM w r 



(2.4) 



In principle, the amplitude should be expanded about the complex poles in k\. Replacing 
the complex pole position by — iM w r w with on-shell mass and width introduces 
an error of 0(a 2 ) (see App. D of Ref. [ |30|1 ). The neglect of the width in the numerator 
introduces an error of order Tw/Mw- Since we apply the DPA only to the corrections, this 
leads to an error of 0(aT^/M^) which is the order of the uncertainty of our calculation. 
Note that additional infrared (IR) singularities appear in the factorizable corrections when 
taking the on-shell limit in the numerator. 

In lowest order, the matrix element in DPA factorizes into the one for the production 
of the two on-shell W bosons, -M|orn^ W+W 5 the (transverse parts of the) propagators of 
these bosons, and the matrix elements for the decays of these on-shell bosons, M^ OT ^^ 1 ^ 2 



and M 



W" 
Born 



M 



Born, DPA 



+WW^4/ 



E 



M 



e T e - 
Born 



- /V/l Born JVl l 



l Born 



(2.5) 



where we introduced the abbreviations 

K± = k\- M 2 , 



M 



Ml- - iM w r 



W-L w 



(2.6) 
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for the off-shellness of the W bosons. Equation (|2.5|) contains the coherent sum over 
the physical polarizations \± of the W ± bosons. Note that the polarizations have to be 
defined in the same way for the production and the decay matrix elements. 

The definition of the DPA for factorizable diagrams ( |2.4|) implies that the squared 
momenta of the W-boson legs of the production and decay vertex functions, which are 
hidden in the shaded blobs of Figure [|, have to be set to their on-shell values M^. 
The W self-energies, marked by open blobs, have to be expanded about M^, and the 
corresponding residues have to be distributed equally to the production and decay parts. 
Finally, we can express the factorizable doubly-resonant corrections by the product of 
gauge-independent on-shell matrix elements for W-pair production and W decays, and 
the (transverse parts of the) W propagators, 

AA/f e+e "^ WW ^ - V 1 f f A >e+e-->W+W- A/ fW+-+/i/ 2 . >W-—/ 3 /4 

0yX/t virt,fact,DPA ~ 2^ T? T? \ 0JVl -^'Born JVl l 



l virt,fact,DPA / j j£ j£ yjjvi -^'Born J l Born 

A + ,A_ + ~ 

+ M^ W+W ~SM W+ ^MZrn M ~' 



+ MZ£-+™-M£?™5M™-+h*), (2.7) 
where 8M. denote one-loop contributions. Details on the actual calculation of 



M 

virt fact dpa can De f° un d in Section |3[I 



2.1.2 Virtual non-factorizable corrections 

All loop diagrams that are not of the generic form of Figure ^ belong to the virtual 
non-factorizable corrections. Simple power counting [ [Tj| reveals that only diagrams with 
photon exchange in the loop can give rise to doubly-resonant (virtual) non-factorizable 
corrections. The considered set of diagrams is not gauge-invariant even in the DPA. There- 
fore, we have defined the non-factorizable doubly-resonant corrections by subtracting the 
factorizable doubly-resonant corrections from the complete doubly-resonant corrections [ 



T6| . The so-defined non-factorizable corrections do not only receive contributions from 
diagrams in which the photon links the production and decay subprocesses, i.e. manifestly 
non-factorizable diagrams [graphs like (a), (b), (c) of Figure 0, but also from factoriz- 
able diagrams [graphs like (d), (e), (f), (g) of Figure 0. The latter contributions arise 
by subtracting the factorizable contributions of those graphs, which contain the artificial 
IR singularities mentioned above, from the complete graphs in DPA, in which the IR 
singularities are replaced by logarithms of the form lnK±. A representative set of dia- 
grams contributing to the non-factorizable corrections to e + e~ — > WW — > 4/ is shown in 
Figure 

Note that the non-factorizable corrections do not contain a product of two indepen- 
dent W propagators and that these corrections involve logarithms of the form ha.(K±), 
which have to be kept exactly. Therefore, it is not possible to define a residue similar 
to R(M$f, Mw, 6) for the double resonance as done for the factorizable corrections. The 
resonance structure of the non-factorizable corrections has the form of a homogeneous 
polynomial of order two in K + and K- in the denominator. 

The non-factorizable corrections yield a simple correction factor S^lLt dpa to the lowest- 



order cross section. Its explicit form is given in Section 3.2 
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(a) type (mf) (b) type (ff) 




Figure 2: A representative set of diagrams contributing to the virtual non-factorizable 
corrections 
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2.1.3 On- shell projection 

The gauge-invariant definition of the factorizable corrections requires the introduction 
of an associated phase-space point with on-shell W bosons for each phase-space point with 
general off-shell W bosons. In our formulas, we leave the form of the on-shell projection 
open and mark all momenta of the associated on-shell point by carets. In particular, the 
on-shell projections of the momenta k + = k\ + k 2 and k_ = fc 3 + /c 4 of the W bosons obey 

k+ = k 1 + k 2 , k- = k 3 + h, k 2 ± = M^. (2.8) 

The actual form of this on-shell projection is not uniquely determined, but involves some 
freedom. However, this intrinsic ambiguity leads to differences of the order of oTw / (7rMw) 
for different versions of the projection, i.e. the ambiguity remains below the desired level 
of accuracy for the DPA. The explicit on-shell projection used in the numerics is given in 
App. [A]. For the on-shell kinematics we define the Mandelstam variables s and t as usual, 

s = {p++p-) 2 = {k + + k^)\ 

i = (p+ - Kf = (P- - M 2 = M&- |(1 - /Joes*), (2.9) 

where (3 = \J\ — 4M(y / s and 6 are the velocity and the scattering angle of the outgoing 
W bosons, respectively. 

2.2 Real photonic corrections 

A substantial part of the 0(a) corrections to four-fermion production is due to real 
photon emission. Since detectors are not able to detect photons that are soft or collinear 
to fermions (other than muons), we have events with visible and invisible photons. The 
former are events in which a photon of finite energy is separated from all other particles in 
the detector, the latter correspond to events in which no separate photon is seen. Events 
with visible photons can be taken into account via the processes e + e~ — > 4/7, for invisible 
photons a careful treatment of soft and collinear singularities is required in addition. 

2.2.1 The bremsstrahlung process e + e~ — > 4/7 

For the bremsstrahlung process e + e~ — > 4/7 the complete lowest-order matrix ele- 
ments Ai e+e ^ 4 ^ 7 are taken into account. More precisely, we employ the matrix elements 
for the 4/ + 7 final states of the CC11 class, which is the minimal gauge-invariant subset 
of graphs containing all doubly-resonant diagrams corresponding to W-pair production. 
We use the explicit results presented in Ref. [ ^8| . 

2.2.2 Invisible photons 

The emission of invisible photons contributes to the 0(a) corrections to e + e~ — > Af. 
They include soft and collinear singularities, which are regularized by an infinitesimal 
photon mass A and small fermion masses mf, respectively. Since the singularities of real 
and virtual corrections are related to each other, these corrections have to be combined 
carefully, in order to avoid a mismatch in the singularity structure. We use two different 



8 



procedures of treating soft and collinear photon emission: one is based on a subtraction 
method, the other on phase-space slicing. The precise implementation of these procedures 
is discussed in Section f|. 

2.3 The master formula 

The total cross section is composed of the following contributions 

J da = J da£ f n ^ f + I dalif^ f + J da^~^\ (2.10) 

Here da^ n ~* 4 ^ is the full lowest-order cross section to e + e~ — * 4/, 

d4^ 4/ = ^d$ 4/ |A<^| 2 (2.11) 

with the corresponding matrix element M.^^^^ an d the four-particle phase-space ele- 
ment d$4j. Similarly, d<r e+e ^ 4 ^ 7 , which describes the real corrections, is the full lowest- 
order cross section to e + e~ — ► 4/7, and da^ t denotes the virtual one-loop corrections. 

Both the virtual and the real corrections involve soft and collinear singularities. These 
singularities are extracted by separating the cross sections into finite and singular parts: 

, e +e--+4/ _ i e+e-->4/ , , e + e --4/ 
UL, virt ~~ uc virt, finite ' uu virt,sing ) 

da e +e -^4/ 7 = da ^fr + da^l^- (2.12) 

The singular part of the virtual corrections factorizes into the lowest-order cross section 
and a simple correction factor. Moreover, the singularities in the real part, do"g in g ^ 4 ^ 7 ; 
can be split off, and the five-particle phase-space element d$4/ 7 can be decomposed into 
d$4/ and the phase-space element of the photon d$ 7 . The integration over d$ 7 can be 
partially performed, and the result can be written as a convolution of a structure function 
with the lowest-order cross section da^^^ . When adding both contributions, all soft 
singularities and all collinear singularities associated with the final state cancel, and the 
resulting "singular" cross section, 

, e+e-^4/ = , e +e--v4/ , , e + e -^4/ 7 

virt+real,sing virt,sing 1 ^^sing > \ 

contains, apart from finite terms, only collinear singularities associated with the initial 
state, i.e. leading logarithms of the form ln(s/mg), at least for inclusive enough observ- 
ables. 

In RacoonWW, the DPA is only applied to the finite part of the virtual corrections, 

e+e-^4/ j c+e-^WW^4/ _ , e+ e -^WW->4/ , e+e"^WW^4/ (0-\A\ 
a<J virt,finite ~~ > a<T virt,finite,DPA — a<J virt,DPA a<J virt,sing,DPA i l^- 14 J 



with da^^si^SpA* 4 ^ defined in Section |Q[ The doubly-resonant virtual corrections are 
composed as follows, 

dff virt,DPA = Ys V [ 2R£ l( jM Born,DPA ) ^virt.fect.DPA J 

4. I /u c+e "^ ww ^i 2 A virt 1 c? 

+ K ' Born, DPA I °nfact,DPA i \ z - LO ) 
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where •A4^ ir ^f a ^.^^* 4 ^ denotes the matrix element for the factorizable virtual corrections 
(Section |3J]) and S^ ct DPA is the factor describing the non-factorizable virtual corrections 
(Section |3~^ ). 

Finally, we arrive at the master formula for the cross section: 

( A — f A e+e ~^ 4 / I f j e+e"^WW^4/ . f j e+e"^4/ / , e+e"^4/7 fr> i R \ 

a<r Bora + / oi(7 virtjfinitCjD p A + / acr virt+real)Sing -+- / acr finite . ^z.ioj 

Since the contribution da^^+reai'fsing * s n °t treated in DPA, the leading-logarithmic pho- 
tonic corrections resulting from initial-state radiation (ISR) are treated exactly in our 
approach. 

The separation of the singularities can be done in different ways. In RacoonWW, we 
have implemented two possibilities to extract the singularities from the real corrections, 



one is based on a subtraction method and discussed in Section 4.1, the other is based 



on the phase-space-slicing method and described in Section L2. Both procedures are 
equivalent. For the extraction of the singularities from the virtual corrections, also two 
choices have been implemented. One is motivated by the subtraction method, the other is 
inspired by the definition of a U(l)-gauge-invariant virtual photon part following Ref. [ 
(cf. Section |4.3[ ). Since the definition of a finite part of the virtual corrections differs by 
finite terms between both cases, and since the singular parts are treated exactly but 
the finite parts in DPA, both approaches differ by terms of the order of the intrinsic 
uncertainty of the DPA. These differences are due to finite terms that are either taken 
into account in DPA or exactly. 

3 Virtual corrections in double-pole approximation 

3.1 Factorizable corrections 

We have calculated the contribution of the virtual corrections to the transition matrix 
elements in DPA, 5M. c vi ^ t ^dpA) m two different ways, the results of which are in 
perfect numerical agreement. In a first approach we made use of the results of Refs. [ |l(| 
and [ [13|] for the pair production and decay of on-shell W bosons, respectively. The 
virtual corrections in DPA are given by (|2.7|) in this approach. When using the existing 
on-shell results, particular care is needed concerning the polarization of the intermediate 
W bosons. In Refs. [ [K| [13| the polarization vectors for the W bosons were introduced 
in a way convenient for each process. In order to keep the spin correlation between 
production and decay subprocesses, it is, however, necessary to take one and the same 
polarization vectors for the production and the decay of a W boson and to perform the 
summation over the W polarizations coherently, i.e. at the level of the unsquared matrix 
element. Hence, the results for on-shell W bosons cannot be taken as black boxes, but 
have to be carefully combined. 

In a second approach we have performed a completely new and independent calcula- 
tion of all relevant one-loop diagrams for the pair production and the decay of on-shell 
W bosons. The Feynman graphs have been generated with FeynArts [ |32| . The actual 
calculation has been done twice using our own Mathematica [ |33| routines, once in the 
't Hooft-Feynman gauge using the Feynman rules and the renormalization scheme of 
Ref. [ ^] and once in the background-field formalism using the results of Ref. [ |55|. The 
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evaluation of the one-loop diagrams follows the methods described in Ref. [ EM in both 



calculations, and the two results are in perfect numerical agreement. In the following we 
briefly sketch the strategy of these calculations which exploits the factorization property 
of the amplitudes. 

3.1.1 Structure of the matrix elements 

The DPA matrix elements are decomposed into a set of so-called standard matrix 
elements (SMEs) M a w which contain the spin structure of the external particles, and 
invariant coefficients F£. Since the relevant diagrams factorize into W-pair production 
and W decay, the SMEs for W-pair production can be chosen as (a = ±) 

f -J*+)£.w ff u(p_), 

- $_)uVf/(p_)(4£i), 



+ - jLK«(p_) \(e*Jj)(eU>-) + (e> + )(e!i+)l (3.1) 



Ml 


= v(p+] 


m 


Ml 


= v(p+] 


Uh 


Ml 


= v(p+) 




Ml 


= v(P+) 




Ml 


= v{p+] 


Uh 


Ml 


= v(p+] 


Uh 


M a 7 


= v{p+, 


Uh 



with "effective W-polarization vectors" 



V2s w A 4 



= -JL- ±- u(k 3 )Yoo-v(h). (3.2) 

The 14 SMEs given in Q3.1p are exactly those defined in Ref. [ pi] , where v(p+) and u{pJ) 
are the Dirac spinors of the massless e + e~ initial state. The chirality projectors are given 
by uj± — |(1 ± 75). The effective W-polarization vectors of ( ^.2| ) are a formal shorthand 
for the W propagators and the tree-level decay matrix elements. Of course, the off- 
shellness of the W bosons has to be kept in these W propagators, where k+ = k\ + ki and 
A;_ = k 3 + fc 4 are inserted, but the W decays are calculated with the on-shell kinematics, 
i.e. the Dirac spinors u(ki) and v(ki) depend on the on-shell-projected momenta. Owing 
to the consistent neglect of fermion masses, the polarization vectors are transverse, 

fc±4 = 0. (3.3) 

This transversality and k± = ensure that the algebraic reduction of the 0(a) cor- 
rections to the W-pair production subprocess to SMEs and invariant coefficients can be 
performed in the same way as for pure on-shell pair production, i.e. without knowing any 
details of the W decays. Explicitly, the decomposition of DPA matrix elements reads 

M%^]**{p+,p-XXA>#-) = jlK{s^Ml{p + ,p-XXX^-)- (3-4) 

n=l 
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Figure 3: Doubly-resonant lowest-order diagrams for e + e — > WW — > 4/ 

It is important to realize that the invariant coefficients F°(s, i) depend only on the scalar 
products of the momenta of the on-shell-projected W bosons, but not on the momenta of 
their decay products. The explicit expressions of the SMEs are listed in App. [B]. 

A few remarks on the choice of SMEs are appropriate. The choice ( |3.1| ) with ( |3.2| ) is 



obtained by taking all products of SMEs for the subprocesses e + e~ — > W + W~, W + 
/1/2, and W" — > fafz for massless fermions (cf. (5.11) and (5.14) of Ref. [ |3H|]), and by 
subsequently replacing £± „£±,*» by the W propagators. Owing to the purely left-handed 
coupling of the W boson and the use of massless fermions, there is only one SME for 
each decay subprocess. Moreover, we have included the lowest-order coupling factor for 
the decays so that the invariant functions F° for the lowest-order amplitudes and the 
loop corrections for the production process are normalized in the same way as for on- 
shell W-pair production. In particular, the non-vanishing coefficients F° Born (s,t) of the 
lowest-order diagrams (see Figure |3|) in DPA are given by 



-^l,Born( S > V 



e 2 



2sir 



Fg^s, t) = -F£ Bom (s, t) = ^ - (l - . (3.5) 

The chosen set of 14 SMEs is not minimal; more precisely, only 12 SMEs are independent. 
The reduction to independent SMEs can be performed by making use of two relations 
between several SMEs that follow from the four-dimensionality of space time (cf. (11.8) 



of Ref. [0). 

The contributions of the virtual corrections to the W-pair production subprocess to 
the invariant coefficients F° can be directly read of from the results of Ref. [ |l(| . The 
contributions of the W decays are obtained by multiplying the corrections to the decay 



matrix elements, which can be extracted for instance from Ref. [ pl| , with the lowest-order 
coefficients F° Born (s, t) of (|3.5|). Note that care has to be taken that the relevant phases 
are taken into account properly. 



3.1.2 Evaluation of coefficient functions 

Finally, we add some remarks on the actual numerical evaluation of the invariant 
coefficients F°(s,t). The expressions for the one-loop contributions to these coefficients 
are rather involved (see Refs. [ [10|, |i~lf) so that the corresponding computer codes are 
lengthy and slow. Moreover, the employed Passarino-Veltman reduction [ [JBJ of tensor 
integrals to scalar integrals breaks down at the boundary of phase space (cos 9 = ±1), 
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necessarily leading to numerical instabilities in the very forward and backward directions 
of the W-production angle 9. In order to solve these problems of CPU time and numerical 
instabilities, we have made use of the fact that all invariant coefficients F£(s,t) depend 
only on 9 for a fixed scattering energy and fixed input parameters, such as masses and 
couplings. Before starting the multi-dimensional Monte Carlo integration over the four- 
particle phase space, we calculate a set of generalized Fourier coefficients 

2/4-1 r+l 

= \_ x dcos9tF:(s,t)P l (cos9) (3.6) 

by performing the integrals with the Legendre polynomials 

1 d' 



x 2 -l) 1 , 1 = 0,1,..., (3.7) 



numerically. Using, for instance, Gaussian integration for these integrations yields suffi- 
ciently precise results without the need to enter the region of numerical instabilities in the 
coefficients F£(s,t). During the Monte Carlo integration, the coefficients are numerically 
reconstructed by making use of the generalized Fourier series 

OO 1 

i£(a,t) = £T<i(*)A(«**)> (3.8) 
i-o 1 

which involves only trivial algebra. Thus, its evaluation is extremely fast. In practice, 
relatively few generalized Fourier coefficients c^ z (s) are needed; for example, taking I = 
0, . . . , 20 reproduces the full calculation of F£(s, i) within roughly six digits for moderate 
scattering angles and LEP2 energies. Note also that the generalized Fourier series remains 
stable in the forward and backward directions where the original evaluation of F£(s,t) 
breaks down. The explicit factor t in flo\6| ) and ( |3.8j ) was included in order to account for 
a i-channel pole in some of the F£(s,t)\ without this factor the expansion is less efficient 



and requires much more terms in the sum over I in ( |3.8j ). 

3.2 Non-factorizable corrections 

In Refs. [ |37|, [16| the non-factorizable corrections have been discussed in detail. 



In these references, however, the virtual non-factorizable corrections have been combined 
with their counterparts involving real photon emission, in order to obtain an IR-safe 
correction that can be discussed separately. In our present approach we do not separate 
the non-factorizable parts from the full real corrections so that we need the virtual non- 
factorizable corrections separately. Here also those graphs must be taken into account 
that are cancelled by their real counterparts [graphs (c), (e), (f), and (g) of Figure^] 
in the approach of Refs. [ [TJ|, [TJ], |37|, |T^| . The complete expression for the virtual non- 



factorizable corrections is given belowQ, closely following the notation and conventions of 
Ref. [|TJ]. 

The virtual non-factorizable corrections can be written in terms of a correction factor 



^Sact dpa to the lowest-order cross section. Analogously to Ref. [ [T6fl , 6™act dpa * s decom- 



posed into contributions that are associated with a pair of final-state fermions 

a 



St,DPA = E E (-l) a+6+1 Q a Q b - Re{A virt (p + ,P-; k + , k a , fc_, k b )}, (3.9) 



=1,2 b=3,4 



4 The complete expression for the virtual non-factorizable corrections can also be found in Ref. [ p 
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and we only give A virt = A virt (p + , p_; k + , k 2 ; k_, k 3 ), since the other terms follow by 
obvious substitutions. The function A virt receives contributions from diagrams that have 
been classified into seven different types in Figure 0, 



Avirt a virt I A virt , a virt , a virt i a virt i a virt i a virt (n i r\\ 

— A mf + A ff + A if + A mm' + A im + A mf + ^mnr (3.1UJ 

We note that in contrast to the sum of the virtual and real contributions given in 



Refs. [ [14], [15], |37|, |16[ the pure virtual non-factorizable corrections depend both on the 
final and initial state of the reaction. The contributions of the different types of diagrams 
are given by 



Avirt 
mf 



{s 23 + s M )K + D (-k„, k+, k 2 , 0, M, M, m 2 ) 
(sia + s 23 )K_D (-k 3 , -k-, k+, 0, m 3 , M, M), 



A^ rt s 23 K + K„E (-k 3 , -k-, k+, k 2 , A, m 3 , M, M, m 2 ) 



(3.11) 
(3.12) 



a virt 



t +2 K + D (p + , k+, k 2 , A, m e , M, m 2 ) + u_ 2 K + D (p_, k + , k 2 , A, m e , M, m 2 ) 
t- 3 K_D (p^, k-, k 3 , A, m e , M, m 3 ) + u +3 K^D (p + , k 3 , A, m e , M, m 3 ) 



A 



virt 



(2M^ - s)|C (fc+, -k-, 0, M, M) - [C (k+, —k-, A, M w , M w ) 

-(t-M^)| (7 (p+,&+,0,m e ,M) - C (p + ,A; + , A,m e ,M w ) 
+ C (p_, k-, 0, m e , M) - C (p-, k-, A, m e , M w ) 



+ (it - M 



W/ 



C (p + , 0, m e , M) - C (p+, fc_, A, m e , M w ) 
+ C (P-, fc+, 0, m e , M) - C (p-, k + , A, m e , M w ) 



(3.13) 
, (3.14) 



,(3.15) 



a virt 
A mf 



C (k + ,k 2 ,0,M,m 2 ) - [C {k + , k 2 , A, M w ,m 2 ) 
+ C (A;-, k 3 , 0, M, m 3 ) - \C {k-, k 3 , A, M w , m 3 ) 



(3.16) 



A 



virt 



2 /5o(^0 1 M)- J Bo(M 2 ,0,M) ( 5 (fc2, 0, M) - ^(M 2 , 0, M) 
^ iW w1 T, 



M 2 



k 2 _-M 2 



-2B' (M* V ,\,M W ) 
with the kinematical invariants 



(3.17) 





— {k>i + kj) 2 , 




t+i 


= (p+ - h) 2 , 


U-i 


t-i 


= (p- - h)\ 




t 


= (P+ - k+) 2 , 


u 



2 

I) i 



(p_ — k 
(P+ 
(P+ " fc-) 2 



ki) 2 , 



i,3 = 1,2,3,4, 



1,2, 
3,4, 



(3.1* 
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The sign "~" in (|3.11 )-( ^.17|) indicates that the limits k± — > and T w — > are 
implicitly understood whenever possible. The definition of the scalar integrals B , Co, 
D , E and of their arguments can be found in Refs. [ [33], The explicit expressions of 
these functions have been given in Ref. [ [L6[ for the mf', ff', and mm' parts; the ones of 
the remaining scalar integrals are listed in App. y. 

In order to facilitate the evaluation of A vlrt as much as possible, we insert the explicit 
expressions for the scalar integrals into the different contributions. For A^, we can 
simply take over the result of Ref. [ specifically, the combination of Cq functions in 
/ is given there in (C.l) and (C.2). For we obtain 



A virt 
mm 



A virt ~ 2 In 

mm 



AM 



w 



2 In 



AM 



w 



-K. 



+ 4 



(3.19) 



using (|C.1|) of the appendix below. For the remaining parts it is convenient to add up all 
contributions, resulting in 



Avirt I a virt , a virt , A virt , a virt 
mf / + A ff , + A if + A im + A mf 



K + K_s 23 det(Fo) 

det(y) 
K + det(Y 3 



D (-k 4 , k + + k 3 , k 2 + k 3 , 0, M, M, 0) 



+ 2 In 



det(y) 



AM 



^3 
111 



k_ det(r 2 ; 



+ Li 2 1 + 



w 



det(Y) 
u_ 2 (M&-t) 



F 2 + ln 



A 2 



In 



t +2 (M^-u) 



t 



-Li 2 1 + 



t+: 



Li 2 1 + 



+ 2 In 



-K. 



AM 



w 



u 



W-2 



U+3 



s 23 



hi 



ie 



u +3 (M^ - t) 
£_ 3 (M^ - u) 



(3.20) 



where the determinants det(F), det(Yi) and the functions D (. . .), are given in Ref. [ |1| 
[see (3.36), (4.10), and (C.3) there]. Finally, we note that in S^ ct DPA all fermion-mass sin- 
gularities of the initial- and final-state fermions drop out, although some non-factorizable 
graphs contain such fermion-mass singularities. 



4 Treatment of soft and collinear photon emission 

In the following we describe two procedures of treating soft and collinear photon 
emission: one is based on a subtraction method, the other on phase-space slicing. In both 
cases soft and collinear singularities are regularized by an infinitesimal photon mass and 
small fermion masses, respectively. 

For convenience, we introduce a second, generic notation for the momenta and helicities 
of the external particles: 

e + (<?i, «i) + e~(g 2 , k 2 ) -> fi(q 3 , k 3 ) + Mq^ K i) + /a (ft, «s) + h(lG, «e) ( + l(k, A 7 )) 

(4.1) 

in addition to (|2.1|) and ( |2.2| ). Moreover, the masses of the external fermions are denoted 
by rrii (qf = mf — > 0), and we use the invariants s^- = in the following. 
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4.1 The dipole subtraction approach 



The idea of so-called subtraction methods is to subtract a simple auxiliary function 
from the singular integrand of the bremsstrahlung integral and to add this contribution 
back again after partial analytic integration. This auxiliary function has to be chosen in 
such a way that it cancels all singularities of the original integrand so that the phase- 
space integration of the difference can be performed numerically, even over the singular 
regions of the original integrand, which is |.M e+e ^ 4 ^ 7 | 2 in our case. In this difference 
A4 C+C ~^ 4 f~t can be evaluated without regulators for soft or collinear singularities, i.e. we 
can make use of the results of Ref. [ for Ai e+e ^ 4 -^ 7 with massless fermions. The 



auxiliary function has to be simple enough so that it can be integrated over the singular 
regions analytically, when the subtracted contribution is added again. This part contains 
the singular contributions and requires regulators, i.e. photon and fermion masses have to 
be reintroduced there. In RacoonWW we have applied the dipole subtraction formalism, 
which is a process-independent approach that was first proposed [ [38| within QCD for 
massless unpolarized partons and subsequently generalized to photon radiation of massive 
polarized fermions [ |27| . We only need the limit of small fermion masses [ |27|, [L7| in which 



the application of the method is relatively simple. In order to keep the description of the 
method transparent, we describe only the basic structure of the individual terms in ( |2.16|) 
explicitly and refer to Refs. [ |27|, [T7|[] for the details. 

In the dipole subtraction formalism the subtraction function is constructed from con- 
tributions that are labelled by ordered pairs ij of charged fermions, so-called "dipoles". 
The fermions % and j are called emitter and spectator, respectively, since by construction 
only the kinematics of the emitter i leads to collinear singularities. The auxiliary func- 
tions |.M su b,r/| 2 that are subtracted from the original integrand |.M e+e ^ 4 ^ 7 | 2 are given 
by 

|-M SU b^($4/ 7 )| 2 = -(-l) i+j QiQ je 2 £ g^fa^k) IMt^ihuj,™^ 2 , (4.2) 

r=± 

where TKi denotes the helicity of the emitter^. Here and in the following the indices i, j run 
over 1, . . . , 6 if not stated otherwise. Note that |-M S ub,r/| 2 is a function on the entire 4/7 
phase space $4/ 7 , while the lowest-order matrix element -M^n - * 4 ^ on the r.h.s. requires 
momenta of the non-radiative 4/ phase space $4/, i.e. $4/ has to be embedded in <3> 4 j- 7 
by an appropriate mapping. This mapping is chosen differently for different ij pairs, as 
indicated by $4/,^ in the argument of -M-Born^^ ■ The mappings have to ensure that the 
kinematics in $ 4 / 7 and $4/,^ asymptotically approach each other in the soft limit, 



fc^O 



and in the collinear limits 



Qi-k k ~ qi, fc^fc, » = 1,2, (4.4) 
K + k k ~ qi, Qj&ZoQj, i = 3,..., 6, (4.5) 



5 The subtraction functions of Ref. [ |21J and Ref. [ [TtJ are slightly different. As default, the subtraction 
function of Ref. [ [27|] are implemented in RacoonWW. 

6 We suppress the polarization arguments as far as possible. 
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as required by the factorization theorems for mass singularities. The momenta qj have 
to respect momentum conservation and mass-shell conditions everywhere. The process- 
independent radiator functions glj^± behave in the soft limit as 



and in the collinear limits as 



(sub) , 



: ^-,r-J qi k~>0 q . k 
(sub), (1 



x. 



1 + x\ m{ 
Xi 2,qik 



Xi) 2 m 2 



Xi 2(qik) 2 ' 



X A 



1 - 



k° 



gi k-*o q . k 

'I - zA 2 



1 + zf mj 
Zi 2,qik 



mi 



Zi 2(q t ky 
where Pff(y) is the usual splitting function, 



q°i+k^ 



p ( \ 1+ y 2 



y 



1,2, 



(4.6) 



(4.7) 



3,..., 6, (4.8) 



(4.9) 



Taking into account charge conservation, it is easy to show that the sum of the subtraction 
functions \M. su b,ij\ 2 and J2x \M C+C ^ 4 ^ 7 | 2 are asymptotically equal in the singular limits, 



£|A4 e+e "^($ 4 / 7 )| 2 ~ E |AW($4/ 7 )| 2 ior k^O or Ql k^0, (4.10) 



so that the integral of their difference becomes integrable in those regions and the photon 
and fermion masses can be neglected. 

Now we are able to define the finite part of the real corrections, 



/ 



d^r 4/7 = 1/ (1<I», , 



EI-M c+c ^ 4/7 l 2 e($ 4/7 )- £ iM^e^ij) 



(4.11) 



where the photon recombination procedure and separation cuts are included in the ob- 
servable defined by 0. We define 6($4/( 7 )) = 1 if the event passes the separation cuts 
after eventual photon recombination and ©(^4/(7)) = otherwise. In contrast to the first 
term, the observable in the second term depends on the 4/ phase space, 3> 4 /,ij, and is inde- 
pendent of the photon recombination procedure. The reason for the different arguments is 
that the final-state mass singularities included in the subtraction functions have to match 
the corresponding singularities of the virtual corrections exactly. On the other hand, the 
different arguments of the terms require that all photons have to be combined with the 
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nearest charged fermion in the collinear limits in order to obtain B($4/ 7 
the soft and collinear limits [see 



111 



D and (PI)]. 

The calculation of distributions is similar to the application of cuts, since a histogram 
of a distribution is nothing but a series of cuts. Hence, the histogram routine that gen- 
erates the desired distribution during the Monte Carlo integration has to handle each 
column of the histogram in the same way as a cutted contribution to the integrated cross 
section. Note that in this procedure the original differential cross section and the sub- 
traction functions may contribute to different columns of the histogram for one and the 
same event. The final result for each column is nevertheless finite, because such events 
are in general far away from the singular regions[]. 

In addition to the finite part of the real corrections, which is given in ( 4.11 ), we have 
to determine the singular part of the real corrections, da^ g ^ 4 ^ 7 . Therefore, we have to 
evaluate the phase-space integral of all contributions |-M su b,ij| 2 at least over the singular 
regions for finite mass regulators m, and A. To this end, we split the five-particle phase 
space into the four-particle phase space and the remaining photonic parts: 



d$ 



4/7 



dx / d$ 



d$ 



1W 



(4.12) 



If % and/or j are initial-state particles, the momentum of one of these is reduced by a 
factor x, i.e. q% = xqi or q 2 = xq 2 . This is indicated by the argument x in $ 4 j j J (x). (For 
final-state emitter and spectator d$ 7i jj includes the function 5(1 — x) which fixes x = 1.) 
Then, the singular part of the real corrections can be written as 



Z7F i,j=l T=± 



a 



6 



x J dx j d$4f t ij(x)gj^\sij,x 



(sub) 



(4.13) 



M£ c rn ^ 4f ($ 4Uj (x), ™«) | 2 e(<§4/,ii(z) 



ator functions gf^f 1 , 



where the functions T originate from the photonic phase-space integral over the radi- 



( sub )/s.. T ) 



8tt 2 J d$ 7iij x gffi {qi,qj,k), 



(4.14) 



and Sij = 2qiqj. The momenta q iy q~j result from mapping the momenta q iy qj, k to $4/ 
The functions involve mass singularities from the initial state, 



7 (sub) / ~ 



smg 



1,2, 



(4.15) 



which are associated to the splitting function. These singular terms are exactly the well- 
known universal logarithms of the structure-function approach for leading-logarithmic 
ISR. 



7 At the edges of the histogram columns this can also occur for "singular events". The finiteness of 
such contributions is guaranteed by the suppression of phase space for those events. 
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For both emitter and spectator from the final state = 3, ... ,6) the convolution 
over x in (|4.13| ) is absent, so that the integrand depends on the 4/ kinematics of the 
original CM system in this case. For the other cases, where at least the emitter or the 
spectator are from the initial state, the soft singularities appear at the endpoint x — > 1 in 
the convolution over x. This IR-sensitive endpoint contributions can be separated from 
the convolution as follows, 



1 c 1 c 
dosing 



,4/7 



a 
2^ 



E E (-i) i+j QiQi 



(4.16) 



dx 



d$ 



4f,ij { x ) Qij, T 



(sub) / ~ 



&ij > % ) 



2xs 



d$ 



4/,y 



(sub) / ~ 



6 



— y y (-i)' i+j Q 4 Q 7 — 



d$ 



A* 
1 

2s 

(sub) , 



Born 



♦4/ 



(* 



4/,y I 



A* 



Born 



♦4/ 



(* 



4/,y I 



, TK, 



e (^4/,ij(i 

" ©(^(l) 



^s 



Born 



♦4/ 



4/J 



If both i > 2 and j > 2, ^" b ^ is proportional to 5(1 — x) and thus does not contribute 
to the second and third line of (|4.16| ). The second term within square brackets in ( |4.16|) 
represents the endpoint contribution (with opposite sign) at x — > 1 that is split off. All 
the endpoint contributions, which are contained in the last term of (|4.16|) , obey the 4/ 
kinematics of the original CM system, since $4/ = $4/^ (1). The functions originate 
from the integral of the functions over x, 



G 



'y ' 



(sub) / ~ \ 



(4.17) 



As mentioned above, the photon mass A and the fermion masses have to be kept finite, 
in order to regularize soft and collinear singularities in ( |4.17| ). Among these integrals of 
course only the non-spin-flip parts receive singular contributions, 

1 



r {sub) ( ~ s 



where 



In 



V s 



In 



A* 

Sij 



1 

2' 



+ ln(^ 



+ C, 



2' 



In* 



m, 



+ - In 

2 



\ s 



■1 j 



(4.18) 



(4.19) 



These are exactly the soft and collinear singularities of the virtual 0(a) corrections with 
opposite sign. The constants Cy, which are specific for the subtraction terms of Ref. [ |2"?fl, 
read 



a 



ab 



IT* 

Y + 2 ' 



G 



ak 



6 



Cka 



1 — i Cm 

2 2 



(4.20) 



TT Z 3 

-y + 2 

with a,b = 1,2 and Z = 3, . . . , 6. 

Owing to their kinematical structure, the sum of all endpoint contributions, i.e. the 
last line of (|4.16|) , can be directly combined with the virtual corrections before integrating 
over the 4/ phase space. 



19 



4.2 The phase-space-slicing approach 

The idea of the phase-space-slicing method is to divide the 4/7 phase space into 
singular and non-singular regions, then to evaluate the singular regions analytically and 
to perform an explicit cancellation of the arising soft and collinear singularities against 
their counterparts in the virtual corrections. The finite remainder can be evaluated by 
using the usual Monte Carlo techniques. For the actual implementation of this well- 
known procedure (see e.g. Ref. [ |39fl) we closely follow the approaches of Refs. [ [41 . 
We divide the five-particle phase space into soft and collinear regions by introducing the 
cut-off parameters S s and 8 C , respectively. We decompose the real corrections as 



do e+e ~ 



— d$ 4/7 |.M 



e + e 



+4/7 I 2 



do- soft + d(T co ii + doj 



finite 



+4/7 



(4.21) 



Here do so ft describes the contribution of the soft photons, i.e. of photons with energies 
< 5 s i/i/2 = AE in the CM frame, and dcr co ii describes real photon radiation outside 
the soft-photon region (E 1 > AE) but collinear to a charged fermion. We define the 
collinear region by 1 > cos# 7 / > 1 — 5 C , where 6* 7 j is the angle between the charged 
fermion and the emitted photon in the CM frame. The remaining part, which is free of 
singularities, is denoted by dog ni ^ c ^ 4 ^ 7 . 

In the soft and collinear regions, the squared matrix element |yV( e+e ^ 4 ^ 7 | 2 factorizes 
into the leading-order squared matrix element | Born ^ 4jf | 2 and a soft or collinear factor. 
Also the five-particle phase space factorizes into a four-particle and a soft or collinear 
part, so that the integration over the photon phase space can be performed analytically. 

In the soft-photon region, we apply the soft-photon approximation to |.M e+e ^ 4 ^ 7 | 2 , 
i.e. the photon four-momentum k is omitted everywhere but in the IR-singular propaga- 
tors. Since we neglect k also in the resonant gauge-boson propagators we have to assume 
E y < AE rV- In this region do e+e "^ 4/7 can be written as [ ||, H 



do soft = dc< r v 4/ -^£ J2 i-i) l+J QiQi 
47r i=i j=i. 



6 6 



i=i+l 



d 3 k 



Qi 



Jl K E t E x2 £7 \ k ii H. 



1 e+e-^4/ a 

dcr 



Uorn E E {~^-) l+:1 QiQj [lii + Ijj ~ 

i=l j=i+l 



with the basic integrals 



Z7T 



d 3 k 



mi 



E^ {kqi)(kqj)' 



(4.22) 



(4.23) 



The explicit expression for the integrals Ijj can be found in Refs. [ 42|, [34] . For our purpose 
it is sufficient to keep the fermion masses only as regulators for the collinear singularities 
(Ei 3> mi). In this limit we obtain 



do, 



soft 



do 



^E E 



Born 



2tt 



/2AE\ 

2k hr 




- 2 In 



AEiE, 
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mt 



+ -ln 2 



'4K 2 



+ 



2tt 2 



+ 2 Li 2 1 - 



AE,E,; 



(4.24) 



Now we turn to the collinear singularities. In the collinear region, we consider an 
initial-state (final-state) fermion with momentum being split into a collinear photon 
and a fermion with the resulting momentum q~i after (before) photon radiation, i.e. for 
initial-state radiation 



e(g<) -> J{k = (1 - Xi)q^j + e(qi = a^ft) , 
and for final-state radiation 

/(ft = t(* = - Zi)/Zi) + /(ft), 



i = l,2, 



3, 



6. 



In the asymptotic limit, |.M e+e ^ 4 ^ 7 | 2 factorizes into the leading-order squared matrix 
element |.M e+e ^ 4 ^| 2 and a collinear factor describing collinear initial-state and final- 
state radiation, respectively, as long as 5 C is sufficiently small. In the collinear region also 
the five-particle phase space factorizes into a four-particle phase space and a collinear 
part, so that the cross section for hard photon radiation (E^ > AE) in the collinear 
region 1 > cos 6 1 f > 1 — 5 C reads 



da coll = do^ ial + dof™ 1 



(4.25) 



with E{ 3> mi, #7/ = 0{rrii/Ei) and 



da, 



initial 
coll 



dxi 



In f-»^± 

V mf 2 x, ; , 



(4.26) 



and 



i=3 



do - ' 
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Born 



e--+4/ 



(ft,+ 



+ dehorn (ft, - K i 



rl-AE/Ei 



dZiPff(Zi) 



dzi (1 - 



In 



'4Ef5 c ; 
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+ do^n ^ (ft, ~ | 



1 - In 
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5 2tt 2 ' 
+ 2"- 



(4.27) 



with the splitting function Pff of (4.9). In the case of initial-state radiation, the four- 
particle phase space is generated in the CM system of the hard scattering process after 
the emission of the collinear photon, i.e. in the system with CM energy s = x\^s, and the 
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four-momenta are then Lorentz-boosted back to the laboratory frame. In the case of final- 
state radiation the CM energy before and after collinear photon radiation are the same, 
so that dcrB r n does not depend on Zj. Thus, the integration over z t can be performed 



analytically, as done in the last equation of ( 4.27|) . Note that this procedure implicitly 
assumes that photons within small cones collinear to charged final-state fermions will 
never be separated from those collinear fermions. 

Subtracting the soft and collinear cross sections ( |4.24| ) and (|4.25|) from the cross section 
of the process e + e~ — > 4/7 ( |4.21|) yields the finite cross section da^ te ^ 4 ^ . As usual in 
the phase-space-slicing approach, this subtraction is done in practice by imposing cuts 
on the 4/7 phase space, i.e. a photon-energy cut, E 1 > AE = 5 s y/s/2, and a cut on the 
angles between the photon and charged fermions, — 1 < cos# 7 / < 1 — 8 C , but by using the 



exact matrix elements for e + e — > 4/7 of Ref. [ 28|. The collinear and soft cross sections 
are added, 

dosing ~* 4/7 = dcr soft + dcr co ii, (4.28) 
and combined with the singular part of the virtual cross section da^sing 4 ^ of Section 



as described in Section 2.3 



The full radiative corrections are split into a 4/ part with no or an invisible photon, 
d <irt,finTt^DPA 4/ + da^ + ^ 8illg , and a 4/7 part with a visible photon, dcXfi ni e te ^ 4/7 . Both 
contributions depend on the cut-off parameters S S ,S C . The dependence on these technical 
cuts cancels in the sum when the cut-off parameters are chosen to be small enough so 
that the soft-photon and leading-pole approximations apply. In Figure ^ we display the 
cut-off dependence of the total cross sections to both parts separately and also illustrate 
the cancellation of the cut-off dependence in their sum. 

4.3 Definition of finite virtual corrections 

Next, we turn to the extraction of the singular contributions from the virtual correc- 



tions, i.e. the definition of dcr^^sh^ in (]2.12|). Here we can make use of the well-known 



cancellation for soft and collinear singularities between real and virtual corrections (KLN 
theorem [ fE|). For inclusive photons, the only uncancelled singularities are the collinear 
singularities originating from the continuum part of initial-state radiation (a^ ^ 1). In 
other words, the singularities of the virtual corrections are exactly given by the endpoint 



parts defined in Section 4.1 within the subtraction approach, but with opposite sign. 



Therefore, we can define the singular part of the virtual corrections by 
do-SjSib = d ^o r e n"^ 4/ ^ E E i-iy^QiQi m?) + C( Sij , m)) + C l3 + Cj 

i=l j=i+l 



3% J 

(4.29) 



with C(sij, mf) and Cy defined in ( [4.19| ) and ( |4.20| ) of Section f£T| , respectively. Note that 
we do not include the polarization-dependent contributions of the endpoint parts [the 
terms ±1/2 in (|4.18| )1 in this definition. Since the polarization-dependent terms are non- 
singular, the simple formulae of this section hold also for the polarized case. Subtracting 
( |1.29| ) from the virtual corrections yields the finite virtual cross section in the subtraction- 
method-inspired approach. We note that it is gauge-invariant by construction, since only 
the Born cross section enters without modification. 
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Figure 4: Cut-off dependence of the total cross section in the phase-space-slicing approach 
for e + e~ — > ud/i - ^ at y/s = 200 GeV. The 4/ part (without Born), the 4/7 part, and 
their sum are shown as a function of log 10 5 S (l.h.s.) (with 5 C = 0.0005) and log 10 8 C (r.h.s.) 
(with 5 S = 0.001). 



The definition of the virtual singular cross section is a matter of convention as far as 
finite contributions are concerned, since finite terms can be redistributed between singular 
and finite parts. Another possibility is to define a finite virtual photon contribution 
by following the approach of Ref. [ [31]. In this YFS-inspired approach, the IR-singular 
contributions are extracted with the help of an explicitly U(l)-gauge-invariant current as 
follows, 



d^YFS 47 = d^ + or V 4/ Re<| (-i4vra)£ £ irV^QiQi 

i=l j=i-\-l 



d 4 k 



(k - 20iq 



(k + 26 jqj ) 



k 2 - 26 iqi k + ie k 2 + 2Q^k + ie 



a 



j 



i=l j=i+l 



A 2 + ie 



(4.30) 



C^Sij.m^ + C(sij,mj) + 2 In 



7r 



^ miirij t 



e(9 l 9 J )n 2 + 2 



where ^ = +1(— 1) if i is a final-state (initial-state) fermion. Again we keep the fermion 
masses only as regulators for the collinear singularities. As can easily be seen from ( [4 .3 Op . 
the U(l) gauge invariance is guaranteed by the existence of a conserved current. This 
YFS factor takes care only of the IR-singular logarithms and the double-logarithmic mass 
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singularities. Since the single-logarithmic mass- singular terms build a separately gauge- 
invariant subset we can simply add them to the YFS factor as follows [ H , 



i e+e -^4/ 
acr virt, sing, YFS 



J e+e" -4/ 
uo soft, YFS 



— da 



Born 



7T 



6 

2. E< 

i=l j'=i+l 



-1 



QiQi In 



(4.31) 



This defines the singular virtual cross section in the YFS-inspired approach. 

While the requirement of the cancellation of the soft and collinear singularities in ( 2.14|) 
unambiguously fixes the logarithmic terms, the two definitions of the virtual singular 
contribution differ by finite, non-logarithmic terms 



j e+e" -44/ 
uw virt, sing, sub 



l e+e"^4/ 
a<T virt,sing,YFS 



(4.32) 
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2vr 



6 6 

E E< 
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7T 



2 + dj + C 



-li- 



as can been seen by comparing (|4.29| ) and (|4.31|) . Since the singular term do"^ irt sing , sub Y fs) 
is treated once in DPA, when subtracted from the doubly-resonant virtual corrections 
d^virt,DRA VW ^ 4 ^ [see ([2.14 )1 , and once exactly, when included in da^ + ^/ si , the cross 
sections calculated with the subtraction-method-inspired and in the YFS-inspired ap- 
proach differ by non-doubly-resonant terms 
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da 



Born 
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, e+e--*WW^4/ 
aCr Born,DPA 



x or E E (-i) i+j QiQj 

Z7F i=l j=i+l 



7T 



+ 9(9 i 9 j )n 2 -2 + dj + C. 



(4.33) 



As can be seen from ( |4.33| ), the ambiguity induced by these finite terms is of the order of 
the uncertainty of the DPA. A numerical discussion of this ambiguity for observables can 
be found in Section [773|. 

While the two different definitions of the virtual singular corrections are inspired by the 
subtraction method and the YFS treatment of IR singularities, their use is independent 
of the method employed for the treatment of singularities in the bremsstrahlung contri- 
bution. This means that for the splitting between singular and finite virtual corrections 
both definitions da^ s ^gf suh Y fs) can ^ e usec ^ * n ^ ne subtraction and phase-space-slicing 
approaches. 



5 Higher-order initial-state radiation 

The emission of photons collinear to the incoming electrons or positrons leads to 
corrections that are enhanced by large logarithms. In order to achieve an accuracy at the 
few 0.1% level, the corresponding higher-order contributions, i.e. contributions beyond 
0(a), must be taken into account. This can be done in the structure-function method [ 
f45| , 0. According to the mass-factorization theorem, the leading-logarithmic (LL) initial- 
state QED corrections can be written as a convolution of the lowest-order cross section 
with structure functions, and the corresponding differential cross section reads 

J da LL = j\ Xl £ dx 2 T^(x 1 ,Q 2 )T^(x 2 ,Q 2 ) J da c ^ 4f (x lP+ ,x 2P ^). (5.1) 
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Here x\ and x 2 denote the fractions of the longitudinal momentum carried by the incom- 
ing electron and positron momenta just before the hard scattering process occurs. This 
means that the incoming momenta p± before emission of the collinear photon are rescaled 
by X\ >2 , and the CM frame of the hard scattering process with the lowest-order cross 
section da^ n (xxp+, x 2 pJ) is boosted along the beam axis. The LL structure function 
including 0(a 3 ) terms is given by [ |7j 
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0e 



2a 



71 



(L-l) 



and the leading logarithm 
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Q 2 



(5.2) 
(5.3) 
(5.4) 



Note that the scale Q 2 is not fixed within LL approximation, but has to be set to a typical 
scale of the underlying process; for the numerics we use Q 2 = s. In (|5.2| ) 7^ is the Euler 
constant and T(y) the gamma function, which should not be confused with the structure 
functions. Note that some non-leading terms are incorporated, taking into account the 
fact that the residue of the soft-photon pole is proportional to L — 1 rather than L for 
the initial-state photon radiation. 

We add the cross section fl5.1| ) to the one-loop result and subtract the lowest-order 
and one-loop contributions da hL ' 1 already contained within this formula, 



da 



LL,1 



dxidx 2 [5(l -x 2 ) + r^ c L ' 1 (xi,Q 2 )5(l - x 2 ) 

+ 5{l-x 1 )T^\x 2 ,Q 2 
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(5.5) 



in order to avoid double counting. The one-loop contribution to the structure function 
reads 
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Note that the uncertainty that is connected with the choice of Q 2 enters now in 0(a 2 ), 
since all 0{ot) corrections, including constant terms, are taken into account. 



6 QCD corrections 

QCD corrections enter the processes e + e~ — > WW — > 4/ in two different places. On 
the one hand, the hadronic W width receives a QCD correction, namely a factor (l+a s /7r) 



33, |13[ m 0{a^). This affects the total W width IV in the resonant W propagators. On 
the other hand, each hadronically decaying W boson receives a QCD correction to the 
Wqq' vertex. If the full phase space for gluon emission is integrated over, this correction 
reduces to a multiplicative correction factor (1 + a s /7r) for each hadronically decaying 
W boson in 0(a s ). The application of this inclusive factor to distributions is usually 
called "naive QCD correction". In the total cross section these naive QCD factors cancel 
against the corresponding factors in the W width in the resonant W propagators. Note 
that one-gluon exchange^] between quarks from different W bosons vanishes exactly for 
pure charged-current (CC) reactions, because of the colour structure. 

If the gluon phase space is not integrated over, QCD corrections have to be calculated 
from the virtual vertex corrections and from real gluon emission with the matrix elements 
e + e~ — > 4/ + g. This option is also supported by RacoonWW for four-fermion final 
states of the CC11 class. The calculation is technically similar to the one for the photonic 
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7 Numerical results 

For the numerical results, with the exception of the comparison to the results of other 



groups, we used (as in Refs. [ ^5], ^6], |J) the following parameters 



G> — 


1.16637 x l(T 5 GeV- 2 , 


a = 


1/137.0359895, 


M w = 


80.35 GeV, 


r w = 


2.08699 . . . GeV, 


M z = 


91.1867 GeV, 


r z = 


2.49471 GeV, 


m t = 


174.17 GeV, 


M H = 


150 GeV, 


m e = 


510.99907 keV. 







(7.1) 

We work in the fixed- width scheme and fix the weak mixing angle by c w = Myj/Mz, 
= 1 — c^. The parameter set ( [7.1Q is over-complete but self-consistent. Instead of a we 
use G M to parametrize the lowest-order matrix element, i.e. we use the effective coupling 

(7.2) 

in the lowest-order matrix element. This parametrization has the advantage that all 
higher-order contributions associated with the running of the electromagnetic coupling 
from zero to and the leading universal two-loop mt-dependent corrections are already 
absorbed in the Born cross section. In the relative 0{a) corrections, on the other hand, 

8 Multiple gluon exchange, which is intrinsically connected to colour reconnection, is not considered 
here. It should be included in the hadronization simulation. 
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we use a = a(0), since in the real corrections and in the mass-singular virtual corrections, 
which yield the bulk of the remaining corrections, the scale of the real or virtual photon is 
zero. The W-boson width given above is calculated including the electroweak and QCD 
one-loop corrections with a s = 0.119. 

As default, we use the following set of separation and recombination cuts: 

1. All photons within a cone of 5 degrees around the beams are treated as invisible, 
i.e. their momenta are disregarded when calculating angles, energies, and invariant 
masses. 

2. Next, the invariant masses Mj 7 of the photon with each of the charged final-state 
fermions are calculated. If the smallest Mj 7 is smaller than a certain cutoff M rec 
or if the energy of the photon is smaller than 1 GeV, the photon is combined with 
the charged final-state fermion that leads to the smallest M/ 7 , i.e. the momenta of 
the photon and the fermion are added and associated with the momentum of the 
fermion, and the photon is discarded^. 

3. Finally, all events are discarded in which one of the charged final-state fermions is 
within a cone of 10 degrees around the beams. No other cuts are applied. 



In Refs. [ Q we have already presented a short survey of numerical results for 
0{a) corrections to e + e" — > WW — > 4/ obtained with RacoonWW for LEP2 and 
linear-collider energies. In the meantime we extended RacoonWW by taking into ac- 
count higher-order ISR and by considering QCD one-loop corrections as described in 
Section [| and Section || respectively. Here we concentrate on LEP2 energies and provide 
results for the total cross sections and distributions obtained with this extended version of 
RacoonWW. Additional numerical results can also be found in Ref. [ [|. Unless stated 
otherwise, for the lowest-order contributions only the CC03 diagrams are taken into ac- 
count. Our "best" results comprise the CC03 Born cross sections and the corrected cross 
sections including all the radiative corrections described in this paper, i.e. electroweak 
one- loop corrections in DPA, exact 0(a) photon radiation based on the full matrix ele- 
ment for e + e~ — > 4/7, higher-order ISR up to 0(a 3 ) and "naive" QCD 0(a s ) corrections 
(see Section ||). 

As mentioned before, RacoonWW involves two branches for the treatment of soft 



and collinear singularities, one following the subtraction (see Section |4.1|) and one the 
phase-space-slicing method (see Section |4~2| ) . While these two branches use the same 
matrix elements, the Monte Carlo integration is performed completely independently, 
thus providing us with a powerful numerical check of RacoonWW. In the following we 
present numerical results of both branches of RacoonWW, starting with the total cross 
sections to e + e~ — > WW — > 4/ at LEP2 CM energies. If not stated otherwise, the finite 
virtual corrections that are treated in DPA are defined by the subtraction-method-inspired 
approach (|4.29|) . 



Except for the 1 GeV cut, the described cut and recombination procedure coincides with the one used 



in Refs. [ 25, E6 
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Figure 5: Total WW production cross section at LEP2, as given by the LEPEWWG [ [IfJ 



7.1 The total W-pair production cross section at LEP2 

In Tables [1] and ^| we list the predictions of RacoonWW for the different W-pair 
production channels as well as for the total CC03 W-pair production cross section for 
LEP2 CM energies. Note that for CC03 and negligible fermion masses the results are 
independent of the final state within these channels. In the calculation of these numbers 
no cuts have been applied. The only difference between the results shown in the two 
tables lies in the ISR beyond 0(a): while Table [I] includes only the 0(a 2 ) contributions 
of ( |5.2| ), Table |2| includes the 0(a 3 ) contributions in addition. There are no significant 
differences between both tables, i.e. the effect of the 0(a 3 ) contributions is smaller than 
the integration error, which is at the level of 0.1%. The given errors are purely statistical. 
The errors for the total cross sections were obtained by adding the (statistically correlated) 
errors of the various channels linearly. 

Figure [| shows a comparison of the results of RacoonWW as presented in Table p] 
(see also Ref. [ ||) and of other Monte Carlo predictions with recent LEP2 data, as given 
by the LEP Electroweak Working Group [ [46| for the Winter 2000 conferences. 

The data are in good agreement with the predictions of RacoonWW and YFSWW3 
At the time of the conferences in Winter 2000 the predictions of YFSWW3 



20, 21 



where about 0.5-0.7% larger than those of RacoonWW, which is somewhat larger than 
the intrinsic DPA ambiguity. Meanwhile, however, the main source of this discrepancy 
was found, and the improved YFSWW3 results differ from the ones of RacoonWW 
only by about 0.3% at LEP2 energies. More details on the conceptual differences of the 
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y/S/ (oreV 


Olcptonic/ 


" semileptonic / J-L) 


Ohadronic/ lI) 


a 1 pb 


172.086 


142.088(71) 


442.50(36) 


1376.14(67) 


12.0934(76) 


176.000 


160.076(78) 


498.03(25) 


1550.04(75) 


13.6171(67) 


182.655 


180.697(89) 


562.22(28) 


1749.48(86) 


15.3708(76) 


188.628 


190.882(96) 


594.31(55) 


1848.07(92) 


16.2420(111) 


191.583 


194.271(118) 


604.12(31) 


1880.19(94) 


16.5187(85) 


195.519 


197.320(123) 


614.11(31) 


1911.45(97) 


16.7910(88) 


199.516 


199.497(103) 


620.53(33) 


1931.28(99) 


16.9670(89) 


201.624 


200.200(104) 


622.65(33) 


1937.94(100) 


17.0254(89) 


210.000 


200.910(107) 


624.95(33) 


1945.00(103) 


17.0876(91) 



Table 1: Racoon WW predictions for W-pair production cross sections for LEP2 CM 
energies including 0(a 2 ) ISR (results for the Winter conferences 2000) 



v/i/ GeV 


Peptonic/ fb 


^semileptonic / ft) 


OWronic/ fb 


a ww /pb 


168.000 


115.639(57) 


359.86(18) 


1120.03(55) 


9.8392(49) 


172.086 


142.043(70) 


442.36(36) 


1375.72(67) 


12.0896(76) 


176.000 


160.192(78) 


498.39(24) 


1551.17(75) 


13.6271(66) 


180.000 


173.511(85) 


539.96(27) 


1679.35(82) 


14.7585(72) 


182.655 


180.669(89) 


562.13(28) 


1749.21(86) 


15.3684(76) 


185.000 


185.474(92) 


576.88(29) 


1794.93(88) 


15.7716(78) 


188.628 


190.959(95) 


594.55(55) 


1848.84(91) 


16.2486(111) 


191.583 


194.271(118) 


604.13(31) 


1880.21(94) 


16.5188(85) 


195.519 


197.435(123) 


614.47(31) 


1912.58(97) 


16.8009(87) 


199.516 


199.640(103) 


620.98(33) 


1932.66(99) 


16.9791(88) 


201.624 


200.271(104) 


622.88(33) 


1938.65(100) 


17.0316(89) 


205.000 


200.742(105) 


624.77(33) 


1943.80(101) 


17.0792(89) 


208.000 


200.973(106) 


625.33(33) 


1945.37(102) 


17.0942(90) 


210.000 


200.888(107) 


624.88(33) 


1944.79(103) 


17.0858(91) 


215.000 


200.334(107) 


623.28(33) 


1938.88(103) 


17.0378(91) 



Table 2: RacoonWW predictions for W-pair production cross sections for LEP2 CM 
energies including 0(a 3 ) ISR (results for the Summer conferences 2000) 
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with bare cuts 


0tot[fb] 


final state 


program 


Born 


best 




slicing 


211.166(36) 


191.847(91) 




subtraction 


211.034(39) 


191.686(46) 




(sub-sli) / sli 


-0.06(3)% 


-0.08(5)% 




slicing 


627.38(11) 


591.11(25) 




subtraction 


627.22(12) 


590.94(14) 




(sub-sli) / sli 


-0.03(3)% 


-0.03(5)% 




slicing 


1864.79(32) 


1821.72(66) 


udsc 


subtraction 


1864.28(35) 


1821.16(43) 




(sub-sli) / sli 


-0.03(3)% 


-0.03(4)% 



Table 3: Total CC03 cross sections when using the subtraction method and the phase- 
space-slicing method at y/s = 200 GeV with bare cuts. The numbers in parentheses are 
statistical errors corresponding to the last digits. 



two generators, as well as a detailed comparison of numerical results, can be found in 
Ref. [ [|. A brief summary is given in Section [7 . 4 . 2| . Figure [5] also includes the prediction 
provided by GENTLE [ 0, which is 2-2.5% larger than those from Racoon WW and 
YFSWW3. This difference is due to the neglect of non-leading, non-universal 0(a) 
corrections in GENTLE. The 1-2% reduction in the W-pair production cross section 
by these effects could already be seen from the results of Refs. [ §|, 0, § , where such 
corrections to on-shell W-pair production were considered. In summary, the comparison 
between predictions of the electroweak Standard Model and LEP2 data reveals empirical 
evidence of non-leading electroweak radiative corrections beyond the level of universal 
effects. 



7.2 Distributions for e + e — > ud^t at LEP2 

Now we turn to results obtained with the cut and recombination procedure described 
at the beginning of Section [7|. We consider the cases of a tight recombination cut M rec = 
5 GeV ("bare") and of a loose recombination cut M rcc = 25 GeV ("calo"). 

In Table |3| we provide the total cross sections for leptonic, semileptonic and hadronic 
processes at the CM energy y/s = 200 GeV with the bare cuts applied. The corrections are 
—9.1%, —5.8%, and —2.3% for the leptonic, semileptonic and hadronic channel, respec- 
tively. The results obtained when using the phase-space-slicing method agree well with 
those of the subtraction method, i.e. within the statistical errors of about 0.05%. For calo 
cuts we find the same good agreement between both branches. The differences between 
the cross sections for bare and calo cuts are below 0.05%; of course, the lowest-order 
results are not affected by the recombination procedure. 

In the following we show a variety of distributions for the semi-leptonic process 
e + e~ — > ud/i~P M at y/s = 200 GeV at lowest order and when taking into account radiative 
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Figure 6: Distributions in the W + (l.h.s.) and W (r.h.s.) invariant masses for a calo 
setup for e + e~ — > ud/i - ^ at y/s = 200 GeV 

corrections ("best"). We also display the relative corrections dcr/d<7Born — 1 (in per cent) 
to illustrate the effect of the radiative corrections. For the relative corrections we always 
provide the results of both branches of RacoonWW, one employing the subtraction and 
one the phase-space-slicing method. Since we use the subtraction-method-inspired defi- 
nition of the finite virtual corrections in both cases, and thus the same DPA definition, 
the results should agree within integration errors. 

The invariant-mass distributions for the W + and W~ bosons are shown in Figure ^| 
for the calo recombination cut. The relative corrections to the W + and W~ invariant- 
mass distributions are shown in Figures ^ and [8], respectively, for both recombination 
cuts. The invariant masses are obtained from the four-momenta of the decay fermions of 
the W bosons after eventual recombination with the photon four-momentum. They are 
particularly sensitive to the treatment of the photons, and thus a strong dependence of 
the relative corrections on the recombination cut can be observed. A detailed discussion 
of the distortion of the invariant- mass distributions of the W bosons and their origin 
has been given in Ref. [ . The distributions obtained from the phase-space-slicing and 
subtraction methods are compatible with each other. The integration errors are larger 
further away from the W resonance because of the decreasing statistics. 

In Figures P^TT] on the l.h.s. we always show the distributions in lowest order and 
including the radiative corrections for the recombination cut M rec = 25 GeV, and on the 
r.h.s. the corresponding relative corrections for both the subtraction and phase-space- 
slicing methods. 

We define all angles in the laboratory system, which is the CM system of the initial 
state. The distributions in the cosine of the W + production angle 6* w + and of the decay an- 
gle 6* w +a, together with the corresponding relative corrections, are shown in Figures |]and 



10, respectively. In both cases the corrections depend strongly on the angle. The distribu- 



tions obtained with the subtraction and the phase-space-slicing branch of RacoonWW 
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Figure 7: Relative corrections to the W + invariant-mass distributions for the bare (l.h.s.) 
and calo (r.h.s.) setup for e + e~ — > ud/i - ^ at y/s = 200 GeV 




75 76 77 78 79 80 81 82 83 84 85 75 76 77 78 79 80 81 82 83 84 85 

M w - [GeV] M w - [GeV] 



Figure 8: Relative corrections to the W invariant-mass distributions for the bare (l.h.s.) 
and calo (r.h.s.) setup at for e + e~ — > ud/i^z/^ at yfs = 200 GeV 
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Figure 9: Distributions in the cosine of the W + production angle with respect to the e 
beam (l.h.s.) and the relative corrections (r.h.s.) for e + e~ — > ud/i - ^ at y 7 ? = 200 GeV 
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Figure 10: Distributions in the cosine of the d decay angle with respect to the W + direction 
(l.h.s.) and the relative corrections (r.h.s.) for e + e~ — > udyU _ P M at y/s = 200 GeV 
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agree with each other within integration errors, i.e. in general to better than 0.5%. The 
increase of the integration errors for large decay angles is due to the smallness of the 
corresponding cross section. 

Finally, we consider distributions in azimuthal angles. In Figure |ll] we show the 
distributions over the azimuthal decay angle W + of the W + boson, i.e. the angle between 
the decay plane of the W + and the plane of W-pair production, 

(k+ x p+)(k+ x ki) . 
cos0w+ = Tj |f; ] — r, sgn(sm0 w +) = sgn {k + ■ [(k+ x p + ) x (k+ x ki)]} . 

| Kj_ X p_|- 1 1 *&-\- X K.^ | 

(7.3) 

In this distribution the corrections vary only weakly with the angle. 

The distributions over the angle between the two planes spanned by the momenta 
of the two fermion pairs in which the W bosons decay, i.e. (note that k + = — k for 
non-photonic events) 

cos0 = ^+ TTTTr ^ ^i^m sgn(sin0) = sgn{k+ • [(k+ x ki) x (-k_ x kg)]} , 



|k + x ki||— k_ x k 



3 



(7.4) 

are presented in Figure The corrections are about —5% except for angles about 
0° or 180°, i.e. if the two decay planes coincide. This can be understood as follows: 
for non-photonic events the angle is related to the azimuthal decay angles 0\y± by 
W + ± W - plus some constant (depending on the precise definition of 0w-)- Since the 
virtual corrections depend on W ± only via sin W ± and cos W ± , which are rather smooth 
functions, the peaks for = 0° and = 180° cannot originate from the virtual corrections. 
If, however, a hard photon is emitted, the finite angle between k + and — k_ leads to a 
finite angle between the two decay planes except for the rare situation where k x and k 3 
are in the plane spanned by k + and — k . Thus, for hard photonic events = 0° and 
<p = 180° are suppressed, and the large negative corrections at these angles are remnants 
of the virtual counterparts of the collinear singularities of the ISR, which are only partially 
compensated by real photon emission. This interpretation is supported by the numerical 
observation that the peaks are mainly contained in the convolution over x in do"^ ir ^ +I ^/ sing . 
If we had used k + instead of — k_ in the definition of 0, these large corrections would be 
absent. 

Theoretically the shown azimuthal-angle distributions are of particular interest, since 
they are sensitive to the imaginary parts of the loop corrections. As can be deduced from 



App. A of Ref. [ [I9|], the contributions of the imaginary parts of the one-loop corrections 
always involve a factor sin0w+ or a factor sin0\y- together with symmetric functions 
in these angles. Consequently, they average to zero if the azimuthal angles of the decay 
fermions are integrated over. This is obviously the case for all quantities discussed so 
far except for the and 0w+ distributions. While we do not find significant effects of 
imaginary loop parts in the distribution, we find a small impact on the 0w+ distribution. 
In Figure |13| we compare the relative corrections to the 0w+ distribution with ( "imag" ) 
and without ( "def " ) imaginary parts. 
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ure 11: Distributions in the azimuthal W + angle (l.h.s.) and the relative corrections 
i.s.) for e + e~ — > ud/i'z/^ at y/s = 200 GeV 
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ure 12: Distributions in the azimuthal <fi angle (l.h.s.) and the relative corrections 
i.s.) for e + e~ — > ud/i~P M at yfs = 200 GeV 
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Figure 13: Relative corrections to the distribution in the azimuthal angle 0w+ with 
("imag") and without ("def") imaginary parts of the loop integrals for e + e~ — > ud/i - ^ 
&ty/s = 200 GeV 



7.3 Discussion of intrinsic ambiguities and of other options 

In order to investigate the accuracy of the DPA quantitatively, we have modified the 
implementation of the DPA and compared the obtained results. The differences should 
be of the order of aT^/My^. Recall that in RacoonWW only the finite (i.e. non-mass- 
singular) virtual corrections are treated in DPA, while real photon emission is based on 
the full e + e~ — > 4/7 matrix element with the exact five-particle phase space. Thus, 
only finite virtual corrections are affected by the following modifications. Specifically, we 
consider four types of options: 

• Different on-shell projections: 



As explained in Section [2.1.3| , one has to specify a projection of the physical mo- 
menta to a set of momenta for on-shell W-pair production and decay[^], in order 
to define a DPA. This can be done in an obvious way by fixing the direction of 
one of the W bosons and of one of the final-state fermions originating from either 
W boson in the CM frame of the incoming e + e~ pair. The default in RacoonWW, 
which is explicitly specified in App. [A], is to fix the directions of the momenta of the 
fermions (not of the anti-fermions) resulting from the W + and W~ decays ("def"). 
A different projection is obtained by fixing the direction of the anti-fermion from 
the W + decay ("proj") instead of the fermion direction. 

Different definitions of finite virtual corrections: 

In Section |] the matching of soft and collinear singularities between virtual and 
real corrections has been explained in detail. Recall that the redistribution of these 



10 This option only illustrates the effect of different on-shell projections in the four-particle phase space; 
if real photonic corrections were treated in DPA the impact of different projections could be larger. 
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singular parts fixes only the universal, singular parts, while the redistribution of 
non-singular parts is a mere convention. Owing to the asymmetric treatment of the 
corrections (finite virtual in DPA, real and mass-singular virtual from full matrix 
elements), different redistributions of non-singular contributions change the result 
by terms of the order (a/7r)(Tw/Mw), which is beyond the accuracy of the DPA. 
As discussed in Section |], RacoonWW allows to choose between two possible 
definitions of the finite part of the virtual photonic corrections to which the DPA 
is applied. As default, the subtraction- method-inspired approach ( |4.29| ) is chosen 



("def") and compared to the YFS-inspired approach ("eik") ( [4.3 1| ). The two treat- 
ments differ in the finite parts of the contribution da^l (see Section [4 . 3|) , which 



virt,sing 

is subtracted from the virtual corrections in DPA and added in its exact form to the 
singular part of the real photon corrections. As can be seen from the difference A 
of ( f4.33|) , the ambiguity originates from terms of the form (a/n) x n 2 x (9(1) which 



are either multiplied with the DPA or with the full off-shell Born cross sections. 
On-shell versus off-shell Coulomb singularity: 

The Coulomb singularity is (up to higher orders) fully contained in the virtual 
0(a) correction. Performing the on-shell projection of the DPA to the full virtual 
correction leads to the on-shell Coulomb singularity, which is a simple factor of 
an/ (2/3). However, since the Coulomb singularity is an important correction in the 



LEP2 energy range and is also known beyond DPA [ |49 |, RacoonWW includes 
these extra off-shell parts of the Coulomb correction as default. This replacement of 
the Coulomb singularity is performed by adding and subtracting the corresponding 
contributions in the virtual non-factorizable corrections, as described in Ref. [ [H| . 
Switching the extra off-shell parts of the Coulomb correction off ("Coul"), yields an 
effect of the order of the uncertainty of the DPA. Note that for CM energies close to 
the W-pair-production threshold, the on-shell Coulomb singularity is not adequate, 
and the difference between on-shell and off-shell Coulomb singularity cannot be 
viewed as a measure of the theoretical uncertainty. 

• Imaginary parts of virtual corrections: 

As default the imaginary parts of the loop integrals are neglected in RacoonWW. 
However, RacoonWW contains an option ("imag") that takes into account the 
imaginary parts. Since all contributions of these imaginary parts are proportional 
to the sine of the azimuthal decay angle of one of the W bosons and otherwise 
involve only symmetric functions of these angles, the imaginary parts drop out in 
distributions where these decay angles are integrated over. The numerical check of 
this cancellations can serve as a consistency check of RacoonWW. Moreover, this 
option might be useful for observables where the imaginary parts do not drop out. 

In the following table and figures the total cross section and various distributions to 
e + e~ — > ud/i~z/ M (7) are compared for the different versions of the DPA defined above in 
the calo setup, i.e. with the loose recombination cut M rec = 25 GeV. We consider the CM 
energies 172 GeV, 200 GeV, and 500 GeV. The used lowest-order cross section is based on 
the full 4/ matrix element. "Naive" QCD correction factors and ISR corrections beyond 
0(a) are not included in the results of this section. 
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clef 


proj 


eik 


Coul 




y/s = 172 GeV 


cx/pb 


400.39(22) 


400.27(22) 


400.03(22) 


403.54(22) 


S/% 





-0.03 


-0.09 


0.79 




y/s = 200 GeV 


cr/pb 


570.10(37) 


569.93(37) 


570.04(37) 


570.85(37) 


S/% 





-0.03 


-0.01 


0.13 




y/s = 500 GeV 


cr/pb 


190.30(20) 


190.28(20) 


190.45(20) 


190.31(20) 


S/% 





-0.01 


0.08 


0.01 



Table 4: RacoonWW predictions for the total cross section to e + e _ — > udyU~z/ Ai (7) for 
the calo setup at different CM energies in various versions of the DPA and the relative 
differences S = cr/(Tdef — 1 



The results for the total cross section are shown in Table f|. We find relative differences 
at the level of 0.1% As expected, the prediction that is based on the on-shell Coulomb 
correction is somewhat higher than the exact off-shell treatment, since off-shell effects 
screen the positive Coulomb singularity. As mentioned above, for the low CM energy of 
172 GeV the difference between on-shell and off-shell Coulomb singularity, which is quite 
large (0.79%), cannot be viewed as a measure of the theoretical uncertainty. While the 
difference between the two on-shell projections is below 0.03%, the effect of including 
different finite terms in the DPA ("eik") is at the level of 0.1% in the considered energy 
range. When switching on the imaginary parts in the corrections we find the same results 
as without imaginary parts ("def"). 

In Figures [T3|, [13, and |ll] we show the differences of the "proj", "eik", "Coul", and 
"imag" modifications to the default version of the DPA for some distributions at y/s = 
200 GeV. For the distribution in the cosines of the W-production angle 6* w + and in the 
W-decay angle #w+d (Figure [M]) the relative differences are of the order of 0.1% for all 
angles, which is of the expected order for the intrinsic DPA uncertainty. For the /z-energy 
distribution, shown on the l.h.s. of Figure [L5J, the differences are typically of the same 
order, as long as is in the range for W-pair production, which is 20.2 GeV < < 
79.8 GeV at yfs = 200 GeV. Outside this region, the four-fermion process is not dominated 
by the W-pair diagrams, and the DPA is not reliable anymore, which is also indicated 
by large intrinsic ambiguities. The r.h.s. of Figure [15] shows the DPA ambiguities for the 
ud invariant-mass distribution. Within a window of 2P\y around the W resonance the 
relative differences between the considered modifications are also at the level of 0.1-0.3%. 
The differences grow with the distance from the resonance point. As expected, we find 
that the contribution of the imaginary parts to the production-angle and decay-angle 
distributions and to the invariant-mass distributions are compatible with zero. 

Finally, we show in Figure ^ the difference between different versions of the DPA 
for the distributions in the azimuthal angle between the decay planes of the two W 
bosons and in the azimuthal decay angle of the W + boson. The ambiguities resulting 
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Figure 14: Differences between different versions of the DPA for distributions in the cosines 
of the production 6* w + (l.h.s.) and decay w +d (r.h.s.) angles for e + e~ — > ud//~ 1/^(7) at 
\fs = 200 GeV, as described in the text 
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Figure 15: Differences between different versions of the DPA for distributions in the /x 
energy (l.h.s.) and in the ud invariant mass (r.h.s.) for e + e~ — > ud/i~z/ At (7) at yfs = 
200 GeV, as described in the text 
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Figure 16: Differences between different versions of the DPA for distributions in the 
azimuthal angles <p (1-h.s.) and 0w+ (r.h.s.) for e + e~ — > ud/z - ^^) at = 200 GeV, as 
described in the text 



from "proj", "eik", and "Coul" are of the same size as for the production- and decay- 
angle distributions considered above. While for the distribution no significant impact 
of the imaginary parts is visible, the effect influences the 0w+ distribution by a relative 
contribution that is roughly given by 0.7% x sin0 w +- 

The energy dependence of the ambiguities in the distributions reflects the one of the 
total cross section, in particular if the deviations are flat as for the pro duct ion- angle distri- 
bution. The differences for the distributions in the W-decay angle # w +a and in the ud in- 
variant mass, which are not flat, are illustrated in Figures |l7]and[l8| at ^/~s = 172 GeV and 
500 GeV, respectively. As already stated above, for the low CM energy of ^fs = 172 GeV 
the difference between on-shell and off-shell Coulomb singularity cannot be viewed as an 
ambiguity. On the other hand, the effect of the Coulomb singularity is completely neg- 
ligible for -^/i = 500 GeV. Note also that the largest ambiguities, which can become of 
the order of 0.5%, appear when the cross section becomes small. The large differences 
between the two different on-shell projections in the # w +d distribution at 500 GeV are 
due to the increasing size and the larger angular variations of the relative corrections at 
higher energies. This uncertainty could be reduced by treating the dominant corrections 
at high energies, which are due to Sudakov-type corrections of the form (a/71) ln 2 (s/M^) 



23 1, not in DPA but in connection with the full e + e~ — > 4/ cross sections. 
The discussed results illustrate that the intrinsic ambiguities of the DPA, as applied 
in RacoonWW, are at the level of a few per mil, whenever resonant W-pair production 
dominates the considered observable. 
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Figure 17: Differences between different versions of the DPA for distributions in the cosine 
of the decay angle w+ i (l.h.s.) and in the ud invariant mass (r.h.s.) for e + e~ — > \xdjji~u^) 
at y/s = 172 GeV, as described in the text 
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Figure 18: Differences between different versions of the DPA for distributions in the 
cosine of the decay angle 6* w +a (l.h.s.) and in the ud invariant mass (r.h.s.) for e + e~ — > 
udfj,~u^(j) at yfs = 500 GeV, as described in the text 
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7.4 Comparison with existing results 

In Ref. [ U we performed a detailed, tuned comparison of RacoonWW with a semi- 
analytical calculation of the complete 0(a) corrections to e + e~ — > WW 41eptons in 
DPA by Beenakker, Berends and Chapovsky [ [Tj| (called BBC in the following) and with 
the Monte Carlo generator YFSWW3 [ [20], [Hj]. Here we restrict ourselves to a brief 
summary of the outcome and refer to Ref. [ |3| for more details. Moreover, we present a 
brief comparison of the real-photonic corrections obtained by Jegerlehner and Kolodziej [ 
|4"8"fl (JK) based on an exact matrix-element calculation. 

7. 4-1 Comparison with Berends, Beenakker and Chapovsky 

The numerical comparison with BBC has been done for the purely leptonic channel 
e + e~ — > v il [L r T~v T . The conceptual differences between RacoonWW and the BBC 
approach are as follows: BBC treat both the virtual and real corrections in DPA, while 
the complete real bremsstrahlung corrections and the universal leading-logarithmic part 
of the 0(a) ISR correction are based on the exact CC11 matrix elements for e + e~ — > 4/7 
and e + e~ — > 4/ in RacoonWW. The virtual corrections are treated in DPA in both 
programs, but RacoonWW employs the full off-shell phase space with an appropriate 
on-shell projection, whereas BBC have chosen the on-shell phase space with decoupled 
invariant masses of the W bosons. These differences are (at least formally) beyond the 
accuracy of the DPA. Thus, if the complete photonic phase space is integrated over and 
virtual corrections are added, both approaches yield results in DPA accuracy. 

For the total cross section, the differences between the two approaches should be of the 
naively expected accuracy of the DPA of 0.5% or better, whenever the DPA is applicable. 
We found that both calculations agree very well above 185 GeV. Below this energy the 
differences in the implementation of the DPA become visible. The main effect originates 
probably from the different treatment of the 0(a) ISR and the phase space. While the 
uncertainty arising from the DPA is enhanced by leading ISR logarithms in the BBC 
approach, this enhancement is absent in RacoonWW. 

For distributions unavoidable differences arise from the definition of the phase-space 
variables in the presence of photon recombination. When defining the momenta of the 
W bosons for angular distributions, BBC choose to assign the photon to one of the 
production/decay subprocesses. The angles are then determined from the resulting W 
momenta and the original fermion momenta. In RacoonWW, all angles are defined 
from the fermion momenta after eventual photon recombination. The two different angle 
definitions lead to a redistribution of events in the angular distributions, which arises, 
in particular, from hard photon emission. For the considered W-production-angle and 
W-decay-angle distributions the differences between BBC and RacoonWW turned out 
to be of the order of 1% at a typical LEP2 energy, with a tendency to increase with 
increasing energy. 

Invariant-mass distributions of the W bosons depend crucially on the treatment of real 
photons, which is done in a fundamentally different way in RacoonWW compared to 
the treatment of BBC. Therefore, it does not make sense to compare these distributions 
between the two programs. Specifically, BBC define the W invariant masses from the 
fermion momenta only, i.e. without photon recombination at all. The resulting shifts in 
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the maxima of the distributions are negative, up to — 77MeV for e in the final state for 



t/s = 184 GeV, and dominated by mass singularities [ |i~9| . On the other hand, studies 
with RacoonWW show that these shifts are in general positive and of the order of 
several lOMeV [ |25], ^6f if photon recombination is taken into account. 

Finally, we compared the photon-energy spectrum E 1 {da/dE 1 ). Since the 
RacoonWW prediction for the E 1 spectrum is not based on a DPA, but on the full 
matrix element for 4/ + 7 production, this comparison can serve as a consistency check 
of the DPA for the real-photonic corrections as used by BBC. We find an agreement be- 
tween the two approaches within ~ 10%, which is of the order of the naive expectation 
for the DPA error of 0(T\y/AE). Although these results illustrate the reliability of the 
naive error estimate of the DPA used by BBC for the real corrections at least for LEP2 
energies, this observation cannot serve as a proof that a DPA for real-photonic corrections 
will work for any observable in any kinematical situation, as e.g. for higher energies. 



7.4.2 Comparison with YFSWW3 

In contrast to the comparison with BBC, where some improvements of the DPA were 
switched off in RacoonWW, in the comparison with YFSWW3 only the input pa- 
rameters, cuts, and recombination procedures are chosen to be the same, but no further 
modifications of the codes are made. Both event generators produce their "best" results, 
which can also be directly compared to the LEP2 data. 

We compared the total CC03 cross sections for leptonic, semileptonic, and hadronic W- 
pair production channels and various distributions for e + e~ — > ud/i - ^ at yfs = 200 GeV 
for the same setup as used in this paper. 

The "best" results for the total cross sections obtained with RacoonWW and YF- 
SWW3 when no cuts, bare or calo cuts have been applied are consistent with each other. 
The relative differences amount to only 0.3%, which is within the expected DPA accuracy 
of 0.5%. 

The comparison of the distributions has been done for both, calo and bare, recombina- 
tion schemes. We found that the W-production-angle and W-invariant-mass distributions 
are statistically compatible with each other, i.e. they agree within 1%. We also compared 
several photon observables, such as the distributions of the photon energy, of the photon 
production angle, and of the angle between the photon and the nearest charged fermion 
and found relative differences at the 10% level. This is not surprising since visible photons 
are treated quite differently in RacoonWW and YFSWW3: while the predictions for 
photon observables by RacoonWW are based on the exact lowest-order matrix element 
for e + e~ — > 4/7, in YFSWW3 multi-photon radiation to W-pair production (within the 
YFS scheme) is combined with 0(a 2 ) LL photon radiation in the W decays (done by 
PHOTOS). 



7.4-3 Comparison with Jegerlehner and Kolodziej 

In Ref. [ JK have evaluated some cross sections and distributions for semi-leptonic 
channels e + e~ — > 4/7 of the CC11 class. The calculation is based on the exact matrix 
element with finite fermion masses, and finite gauge-boson widths are introduced in the 
so-called fixed-width scheme. The presented bremsstrahlung corrections to the total cross 
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y/s/ GeV 


189 


500 


A£ 7 / GeV 


JK 48 


1285(1) 


571(1) 


0.001 




1285.0(8) 


570.0(6) 


0.1 




1286.3(8) 


570.1(5) 


1.0 


RacoonWW - slicing 


1286(3) 


575(3) 


0.001 




1286.2(4) 


570.4(3) 


0.1 


RacoonWW - subtraction 


1284.9(5) 


570.4(3) 





Table 5: Comparison of the total cross section cx real / fb for e + e 
10~ 6 GeV, as obtained by JK [ H§ and RacoonWW 



ud/i with A 



sections of e + e~ — > 4f, which have been regularized with a small photon mass A, can be 
compared with predictions made with RacoonWW. To this end, the virtual corrections 
are omitted in RacoonWW and consequently no DPA is needed. The real correction 
provided by RacoonWW corresponds to the full e + e~ — > 4/7 cross section in the suc- 
cessive asymptotic limits A — * and mf — *■ 0. The relative difference to the results of 
JK should thus be of the order of m//Mw, where m/ is the largest fermion mass of the 
process. Table |5] shows a comparison of results for the process e + e~ — > ud/i _ P Ai 7, where 



the input parameters of Ref. [ 48 are used consistently, and the photon mass is set to 



A = 10 -6 GeV. The parameter AE 1 represents the cutoff that separates the domains of 
hard and soft photon emission in the phase-space slicing approach. The matching of hard 
and soft photons is correct up to terms of 0(AE 1 ), i.e. the exact results are formally 
reproduced for AE^ — > 0. Since the 4/7 matrix elements in RacoonWW are based 
on massless fermions, there is also a small acollinearity separation cut S c on all photon 
emission angles in the slicing approach; in the numerics we have set S c = 0.001. The 
inclusion of the collinearity regions, which involves finite fermion masses, is described in 
Section O. The table shows acceptable agreement of the slicing and subtraction branches 
of RacoonWW with the results of JK. 



8 Conclusions 



A proper treatment of electroweak radiative corrections to the W-pair production pro- 
cess e + e~ — > WW — > 4/ is mandatory in order to account for the experimental accuracy 
at LEP2 and future linear colliders. A practical calculation of the 0(a) corrections to 
the full e + e~ — > 4/ process is beyond present possibilities. However, as long as W-pair 
channels dominate, the combination of full lowest-order predictions for the four-fermion 
processes e + e~ — > 4/ with those corrections that are enhanced by two resonant W bosons 
approximates the full (9(o;)-corrected cross sections within a relative accuracy of 0.5%, at 
least for not too high energies. A strategy for such a double-pole approximation (DPA) 
for the 0(a) corrections and its implementation in the Monte Carlo event generator 
RacoonWW are described in this paper. 

In RacoonWW only the non-leading virtual 0(a) corrections are treated in DPA, 
while real photonic corrections are based on full matrix elements for the radiative processes 
e + e~ — > 4/7. In this way, potential problems in the definition of a DPA for real photon 
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emission are avoided, and possible DPAs for the real corrections can be controlled. We 
have shown explicitly how existing results for the corrections to on-shell W-pair production 
and W decay can be employed in the DPA; they define the class of so-called factorizable 
corrections. For the remaining, non-factorizable contributions, which are due to photon 
exchange between the various production and decay subprocesses, we have completed 
the results already given in the literature. Since virtual and real corrections are not 
treated on equal footing, particular care is needed in the treatment of soft and collinear 
singularities. The proper cancellation of these singularities has been discussed in detail; in 
practice two different methods (a subtraction method and phase-space slicing) have been 
applied. Finally, the leading contributions from initial-state radiation beyond O(a) have 
been included in RacoonWW, rendering this program a state-of-the-art calculation. 

For LEP2 energies at detailed discussion of numerical results is presented, includ- 
ing total cross sections, angular and invariant-mass distributions. The consistency of 
RacoonWW is shown by comparing the results obtained with the subtraction method 
with the ones of the slicing approach. Moreover, the reliability of the applied DPA is 
demonstrated by various modifications of the DPA within its intrinsic freedom. The ob- 
tained intrinsic ambiguity supports the expectation that the theoretical uncertainty of 
the used DPA is of the order of 0.5% for energies between 170 GeV and 500 GeV. Above 
500-1000 GeV, leading-logarithmic electroweak corrections of higher orders have to be 
taken into account. 

As a first phenomenological application of RacoonWW the predictions for the total 
W-pair production at LEP2 have been compared with LEP2 data and older "improved 
Born predictions" of GENTLE, which were previously used in the data analyses. The 
RacoonWW results agree with the results of the new generator YFSWW3 but lie 2- 
2.5% below the GENTLE predictions. The fact that the experimental LEP2 data favour 
predictions of RacoonWW and YFSWW3 can be viewed as empirical evidence for 
non-leading electroweak radiative corrections beyond the level of universal effects. 
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Appendix 

A Explicit form of an on-shell projection for the four-fermion final state 

When performing the projection to on-shell W-boson momenta, care should be taken 
that the projected momenta lie in the physical region. This can be ensured by fixing 
angles while performing the limit k± — > M^. An obvious way is to fix the direction of one 
of the W bosons and of one of the final-state fermions originating from either W boson 
in the CM frame of the incoming e + e~ pair. We fix the direction of the W + boson, of the 
fermion fx, and of the fermion fa. Thus, we find for the on-shell-projected momenta 
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— 



2k+k x 
M 2 

2k.k* 



UP — h.P _ L.P 

^2 ~ "'1 ) 

%P —%P _ If* 



(A.l) 



with (3 



4M^/s. 



B Explicit expressions for standard matrix elements 



The SMEs M°, which have been defined in Section |3.1| , have been calculated in the 

|50fl . Since we 



Weyl-van der Waerden spinor formalism using the conventions of Ref. 
consistently neglect fermion masses, the final-state fermions are polarized according to 
o\ = — o 2 = o"3 = — o"4 = 5 and the index a determines the polarization of the incoming 



e according to a 



<7_ 



-a + . The explicit expressions for the two sets Mi read 



Mf 

Mt 

Mt 

Mi 
Mi 
Mi 



Mi 
Mi 



Mi 



Mi 

M 2 

M^ 

Ml 

Ml 
M 6 

Mr 



-AA(j> + k 2 )*(j>„h){(p + k±y(j>+k x ) + (hhnhh) 

2Aikkyikk)(ip+k)*{P-k) + ( P+ k)*(pM), 
2A(p + ky{p-k){ikk)*ikk) + ikk)*ikk)), 
-2A( P+ k)^P-k)((kky(kk) + (kk)*(kk)), 
2A( P+ ky(p-ky(p-k}(p-k), 

-2A(p + k 2 y(p+ky(p + k)(p-k), 

A((kky(kk) + (kky(kk))((kky(kk) + (kky(kk) 
x (( P+ ky(p-k) + (p+ky(p-k)), 
A( P+ ky(pM*(p+k)(pM((p+ky(p-k) + (p + k)*{p-k)), 
A^ky(p-k){(kk)*(kk) + (kky(kk) 

x (( P+ ky{p-k) + {p+ky(p-k)), 

A(p+ky(p+k)((kky(kk) + (kkyikk) 

x ((p+ky(p-k) + (p+ky(p-k) 



-AA(p_hy(p + ki){{p+k 2 y( P+ k 3 ) - (hhyihh) 

2A(kky(kk)((p-ky(p+k) + ( P -ky(p+k)), 
2A{p„ky{p + k) ((kk)*(kk) + (kkyikk)), 
-2A(pSk,y( P+ k)(ikkyikk) + ikk)*ikk)), 
2A(p„kyiP-kyiP+k)iP-k), 

-2A( P+ kyiP-kyip+k)iP + k), 

A(ikk)*ikk) + ikkyikk))(ikk)*ikk) + ikk)*ikk) 
x Up-kyiP+k) + ip-k)*ip+k)), 
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M 8 = A(p + k 2 )*(pM)*<j>+h)<j>-fa)[(p-ki)*<j>+ki) + (p-k 2 )*(p + k 2/j . 
= Mp-ky(p-k)((kky{kk) + (kkyikk)] 
x (fp-kytp+k) + (p-hTiP+k)), 
M w = A(p + k)*(p + h)((kky(kik) + (kk)*(k 2 k 3 ) / 

x ({ P _ky(p + k) + (j>-ky{p+k)), (b.i) 

where (pfc) denotes the spinor product for two light-like momenta p M and k^, and A is the 
global factor 



.4 



1 



1 



24 k\-M^ + iM w r w E-m^ + iM w r w ' 

C Scalar integrals for the virtual non-factorizable corrections 



(B.2) 



In Section |3.2| we have given the virtual non-factorizable corrections in terms of scalar 
one-loop integrals. In this appendix we list the explicit expressions for those scalar inte- 
grals that are not already given in Ref . [ |16| : 



D {p±, k±, h, A, m e , M, mi 



2 In 



m e rrii 
—t±i 



In 



AM 



w 



In 1 



/ m e M 



w 



mi 



M 



w 
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T -L 12 (, 



t-M%f 



Do(p±, k T , h, A, m e , M, m*) ~ D (p±, k±, k h A, m e , M, mi 



C (p±, fc±, 0, m e , M) - Cb(p±, fe±, A, m e , M w ) 



In 



m e M w 



In 



#4 



+ In 



-#4 



A 2 



+ In 



M 



w 



7T 2 



Co(p±, k T , 0, m e , M) - C (p±, fc T , A, m e , M w ) 

Co(p±j fc ±, 0, m e , M) - [Co(p±, fc±, A, m e , M w ] 



C (k±, k h 0, M, mi) - C (k±, k h A, M w , 



4=< 



In 



2 

Ml 



In 



-#4 



AM 



w 



+ ln' 



M 



w 



+ 



7T 



6 J' 



/j (£: 2 ,0,M)-/j (M 2 ,0,M) 
4 -M 2 



B' (M^,X,M W ) 



In 



AM 



w 



+ 1 . (C.l) 



Recall that all expression are given for k\ -> M^ and T w -> 0. 
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